arXiv:1502.01628vl [astro-ph.EP] 5 Feb 2015 


Identification and Dynamical Properties of Asteroid Families 


David Nesvorny 

Department of Space Studies, Southwest Research Institute, Boulder 

Miroslav Broz 

Institute of Astronomy, Charles University, Prague 


Valerio Carruba 

Department of Mathematics, UNESP, Guaratingueta 

Asteroids formed in a dynamically quiescent disk but their orbits became gravitationally 
stirred enough by Jupiter to lead to high-speed collisions. As a result, many dozen large 
asteroids have been disrupted by impacts over the age of the Solar System, producing groups of 
fragments known as asteroid families. Here we explain how the asteroid families are identified, 
review their current inventory, and discuss how they can be used to get insights into long-term 
dynamics of main belt asteroids. Electronic tables of the membership for 122 notable families 
are reported on the Planetary Data System node. See related chapters in this volume for the 
significance of asteroid families for studies of physics of large scale collisions, collisional 
history of the main belt, source regions of the near-Earth asteroids, meteorites and dust particles, 
and space weathering. 


1. INTRODUCTION 

As witnessed by the heavily cratered surfaces imaged by 
spacecrafts, the chief geophysical process affecting aster¬ 
oids is impacts. On rare occasions, the impact of a large 
projectile can be so energetic that the target asteroid is vio¬ 
lently torn apart, and the pieces are thrown into space. The 
sites of such cosmic accidents are filled with debris that 
gravitationally accumulate into larger conglomerates, and 
drift away at speeds that are roughly commensurate with 
the escape speed from the original target body (T4sc)- Ini¬ 
tially, all orbits are similar, because I4sc Krb. where 
Vorh — 15-20 km s“^ is the orbital speed of main belt aster¬ 
oids. On longer timescales, however, the orbits are altered 
by gravitational perturbations from planets, and the orbital 
elements of individual bodies start to diverge. 

It may therefore seem challenging to identify fragments 
of a catastrophic collision that happened eons ago. Fortu¬ 
nately, starting with the pioneering work of K. Hirayama 
(Hirayama, 1918; see also Cimrman, 1917), astronomers 
have developed various methods to deal with this issue 
(Section 2). Roughly speaking, these methods consist in 
a transformation that brings the orbital elements at the ob¬ 
served epoch to a standard, called the proper elements 
(Knezevic et al., 2002), that is unchanging in time (or, 
at least, would be unchanging if chaotic dynamics, non- 
gravitational forces, and other perturbations could be ig¬ 
nored). Thus, ideally, daughter fragments produced by 
breakup of a parent asteroid will appear as a group in space 
of the proper elements even gigayears after the original col¬ 
lision. These groups are called asteroid families, or dynam¬ 
ical families to emphasize that they have been identified 


from dynamical considerations. 

Telescopic surveys such as the Sloan Digital Sky Survey 
(SDSS), Wide-field Infrared Survey Explorer (WISE) and 
AKARI All-Sky Survey provide a wealth of data on physi¬ 
cal properties of the main belt asteroids (Ivezic et al., 2001; 
Mainzer et al., 2011; Usui et al., 2013). They have been 
used to cross-link the color and albedo measurements with 
the lists of dynamical families, in much the same way the 
spectroscopic and taxonomic data have previously been ap¬ 
plied to this purpose (see Cellino et al., 2002 for a review). 
This work is useful to physically characterize the asteroid 
families (see chapter by Masiero et al. in this volume), in¬ 
cluding cases where two or more dynamical families over¬ 
lap, and identify distant “halo” family members that would 
otherwise be confused with the local background (e.g., Broz 
and Morbidelli, 2013). Given that the SDSS and WISE cat¬ 
alogs now contain data for over 100,000 unique asteroids, 
it has also become practical to conduct search for families 
in extended space, where the color and/or albedo data are 
taken into account simultaneously with the orbital elements 
(e.g., Parker et al., 2008; Masiero et al., 2013; Carruba et 
al., 2013a). 

The physical data can be used to identify interlopers. 
The problem of interlopers arises because the clustering cri¬ 
terion applied to identify the dynamical families is only a 
rough expression of the true membership. Unrelated as¬ 
teroids that just happen to have nearby values of proper 
elements will be grouped together with the true members, 
and will thus appear in the lists of dynamical families ob¬ 
tained from the proper elements (e.g., Migliorini et al., 
1995). These interlopers, especially the large ones, intro- 
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Fig. 1.— (a) Clustering algorithm applied to the asteroid belt separates dynamical families (yellow) from the background (red), (b) 
Variation in reflectance properties of main belt asteroids. Here we plot ~25,000 asteroids that were observed by both SDSS and WISE. 
The color code was chosen to highlight the albedo/color contrast of different families. 


duce ambiguity in the interpretation of impact conditions 
that produced individual families, and in the implications of 
these studies for asteroid interiors (cf. Michel et al., 2003; 
Nesvorny et al., 2006a). Here we discuss how large inter¬ 
lopers and true largest members in families can be found 
by applying the “V-shape” criterion, which is based on the 
notion that large fragments are ejected at low speeds, and 
have less mobility due to the Yarkovsky effect (Section 4). 

Based on the synthesis of asteroid families extracted 
from recent publications (Mothe-Diniz et al., 2005; Nes¬ 
vorny et al., 2005; Gil-Hutton, 2006; Parker et al., 2008; 
Nesvorny, 2010, 2012; Novakovic et al., 2011; Broz et al., 
2013; Masiero et al., 2011, 2013; Carruba et al., 2013a; Mi¬ 
lan! et al., 2014), we attempt to build a consensus that could 
serve as a starting point for future studies. We classify the 
asteroid families into notable cases (those that have a high 
statistical significance, and are thus real, and/or are notable 
for other reasons) and candidate families (less interesting 
cases where the statistical significance is low or cannot cur¬ 
rently be established). The distinction between notable and 
candidate families is somewhat arbitrary, and will intention¬ 
ally be left strictly undefined, because that is the nature of 
things. We expect that many candidate families will be con¬ 
firmed with more data, and that a few notable families may 
fall into oblivion. The lists of notable families are being 
made available at the Planetary Data System (PDS) node, 
and are discussed in Section 7. 

A new and exciting development in the past decade was 
the detection of several asteroid families with very young 
formation ages. For example, the Karin family was shown 
to have formed only 5.8 ± 0.2 m.y. ago (Nesvorny et al.. 


2002a). These cases are important, because various col- 
lisional and dynamical processes had little time to act on 
these families to alter their properties. The young families 
have thus attracted much attention from scientists studying 
impact physics, space weathering, debris disks, etc. As we 
explain in Section 3, the age of a young family can be deter¬ 
mined by numerically integrating the orbits of its members 
backward in time and demonstrating that they converge to 
each other at some specific time in the past. This is the time 
of a breakup, and the family age, tage, is the time elapsed 
from the breakup event. 

The method of backward integration of orbits only works 
for the families with tage ^10 This is because dynam¬ 
ics of main-belt asteroids on longer timescales is governed 
by chaos, encounters with (1) Ceres and other large aster¬ 
oids, and non-gravitational forces. A complementary sta¬ 
tistical method for the estimation of family age has been 
developed in Vokrouhlicky et al. (2006a,b). The method 
tracks, in detail, how the family structure in semimajor axis 
changes over time as the family members drift away by the 
Yarkovsky effect (Section 5). The semimajor axis spread 
of an older family will generally be greater than that of a 
younger family. A compilation of formation ages of the 
asteroid families can be used to constrain how the popu¬ 
lation of the asteroid belt collisionally evolved over time 
(e.g., Bottke et al., 2005a,b; Cibulkova et al., 2014), and 
how asteroid surfaces age by space weathering (chapter by 
Brunetto et al. in this volume). 






2. IDENTIFICATION METHOD 

Here we discuss the standard method to identify aster¬ 
oid families. This method consists of the (1) computa¬ 
tion of proper elements, or other elements unchanging with 
time, for asteroids with well-known orbits, (2) identifica¬ 
tion of concentrations or groups of asteroids in proper ele¬ 
ment space, and (3) establishing the statistical significance 
of identified groups. These steps are discussed in Sections 
2.1, 2.2 and 2.3. In Section 2.4, we examine the “overlap 
problem” where two or more families overlap in proper ele¬ 
ment space, and need to be separated. For this it is useful to 
consider families in extended space with physical data be¬ 
ing included in addition to the proper elements. The search 
in extended space can also lead to the identification of new 
families (Section 2.5). Very young families (fage ^ 1 m.y.), 
for which the member orbits have not had time to differen¬ 
tially precess away from each other, can also be identified 
as groups in space of the osculating orbital elements (Sec¬ 
tion 2.6). 

2.1 PROPER ELEMENTS 

The ejection speeds of sizable fragments produced by 
collisional breakups of main belt asteroids are generally 
much smaller than their orbital speeds. The fragments will 
therefore initially cluster near the original orbit of their par¬ 
ent body, and will appear as such if the subsequent effects 
of planetary perturbations are removed by projecting orbits 
into space of proper elements. The three most useful proper 
elements are: the proper semimajor axis (ap), the proper 
eccentricity (ep), and the proper inclination (ip). They are 
close equivalents of their osculating element counterparts 
in that they define the size, elongation and tilt of orbits (see 
Note 1 in Section 9). 

The definition of proper elements as quasi-integrals of 
asteroid motion, and the methods used to compute them, 
were explained in the Asteroids III book (Knezevic et al., 
2002). As these definitions and methods have not changed 
much, we do not discuss them here in detail. In brief, the 
proper elements are obtained from the instantaneous oscu¬ 
lating orbital elements by removing periodic oscillations 
produced by gravitational perturbations of planets. This can 
be done analytically, using perturbation theory (Milani and 
Knezevic, 1990, 1994), or numerically, by integrating the 
orbits and applying the Fourier analysis (Knezevic and Mi¬ 
lani, 2000; Knezevic et al., 2002) (Note 2). 

The computation of analytic proper elements is rela¬ 
tively CPU inexpensive. They are made publicly available 
by A. Milani for both numbered and unnumbered multi¬ 
opposition asteroids at the AstDyS node (Note 3). The 
analytic proper elements lose precision for highly-inclined 
orbits, because the expansion of the gravitational potential 
used to calculate them has poor convergence for high incli¬ 
nations. The more-precise synthetic proper elements (pre¬ 
cision generally at least 3 times better than that of analytic 
elements), on the other hand, require a much larger CPU 


investment, and are only made available at the AstDys node 
for the numbered asteroids. 

2.2 CLUSTERING ALGORITHM 

When these methods are applied to the asteroid belt, 
things are brought into focus with dozens of obvious 
clumps, asteroid families, emerging from the background 
(Figure [T^). To identify an asteroid family, researchers 
apply a clustering algorithm to the distribution of aster¬ 
oids in (ap, ep, ip) space. The most commonly used algo¬ 
rithm is the Hierarchical Clustering Method (HCM; Zap- 
pala et al., 1990) (Note 4), which defines a cutoff dis¬ 
tance, dent, and requires that the length of the link be¬ 
tween two neighboring orbits clustered by the algorithm 
is d = (i(ap, ep, ip) < dcut- A common definition of dis¬ 
tance is = (nap)^(/Ca((5ap/ap)^+/Ce(^ep)^+/ci(Jip)^), 
where n is the orbital frequency, ((5ap, (5ep, 6 sin ip) is the 
separation vector between orbits in 3D space of proper ele¬ 
ments, and {ka,ke,ki) are coefficients of the order of unity. 

The main advantage of the HCM over other methods 
is that there is no strong assumption built into the HCM 
about the shape of an asteroid family in proper element 
space. This is because the chain created by linking nearby 
orbits can track down family members even if their orbits 
dynamically evolved to produce an unusual overall shape. 
A prime example of this is the case of the Koronis fam¬ 
ily which is split into two parts by the secular resonance 
g -f 2^5 — 3^6 = 0 at ^2.92 AU, where g is the apsidal 
frequency of an asteroid, and g^ and g^ and the 5th and 6th 
apsidal frequencies of the planetary system. The part of the 
Koronis family with ap > 2.92 AU has larger eccentricity 
than the part with ap < 2.92 AU, because family mem¬ 
bers drifting by the Yarkovsky effect from ap < 2.92 AU 
have their eccentricities increased by interacting with the 
g ^ 2gs — Sge = 0 resonance (Bottke et al., 2001). 

The main disadvantage of the standard HCM, which be¬ 
comes increasingly difficult to overcome with inclusion of 
numerous small asteroids in the new catalogs, is the prob¬ 
lem of chaining. This problem arises because small frag¬ 
ments are typically ejected at higher speeds and have larger 
mobility due to the Yarkovsky effect. They therefore spread 
more, tend to be distributed more homogeneously through¬ 
out the main belt, and create bridges between different fam¬ 
ilies if a single (large) value of d^nt is used. Clearly, d^nt 
should be set proportional to the asteroid size, or inversely 
proportional to the absolute magnitude H. Expressing this 
dependence, however, adds additional parameters to the 
HCM and makes the whole identification procedure more 
complex. Therefore, in reality, it is preferred to bypass the 
problem of chaining by artificial means (e.g., cuts in proper 
element space applied to deal with individual cases), or the 
proportionality is approximated by a two-step method with 
different cutoffs for small and large bodies (Milani et al., 
2014). 

A tricky part of the HCM algorithm is the choice of the 
cutoff distance. If the value of dent is too small, many 
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dispersed but real families will remain unnoticed and large 
families will artificially be split into parts. If the value is too 
large, the algorithm will clump different families together, 
and will identify irrelevant clumps produced by random 
fluctuations. While many asteroid families can be identi¬ 
fied for a wide range of dcut values, and are real beyond 
doubt, some cases require a specific choice of dcut and can 
potentially be confused with random fiuctuations. Clearly, 
the statistical significance of the identified groups, or their 
insignificance, needs to be established before proceeding 
further. 

2.3 STATISTICAL SIGNIFICANCE 

To make reasonably sure that families identified from the 
HCM are real one can opt for a conservative choice of dcut- 
This can be done, for example, by collecting asteroids in 
a given region of proper element space and redistributing 
them randomly in that region. The HCM applied to this arti¬ 
ficial distribution will reveal that the largest identified group 
with given dcut contains A^*(dcut) members. Now, this pro¬ 
cedure can be repeated, say, one thousand times, recording 
the largest (dcut) obtained from these trials. We can then 
be 99% confident that random fluctuations cannot produce 
groups with more that (dcut) members (this conservative 
estimate includes a 99% confidence interval computed by 
the Wilson score interval approximation). Any group iden¬ 
tified in the real distribution with N > A^*(dcut) members 
is therefore reasonably likely to be real. Higher confidence 
levels can be achieved by increasing the sample size. 

This basic concept, and various modifications of it, is 
known as the Quasi Random Level (QRL; Zappala et al., 
1994). In the ideal world, the QRL would be the ultimate 
solution to the family identification problem: just choose 
dcut and pick up all clumps with more than N > A^* (dcut); 
those clumps are real. Then, there is the real world. First, 
the number density in proper element space is variable due 
to the primordial sculpting of the main belt and resonances 
(e.g., Minton and Malhotra, 2009). Applying a global QRL 
value in (parts of) the main belt may therefore lead to un¬ 
satisfactory results. Second, families do not live in isola¬ 
tion but are frequently close to each other, overlap, and/or 
are surrounded by empty regions. This introduces an ambi¬ 
guity in the QRL definition, because it is not clear a priori 
what region in (ap, ep,ip) space should be considered to 
define the local QRL in the first place. Results may depend 
on this choice. 

In practice, the first choice made is often the minimum 
number of group members, A/min, that is considered to 
be interesting. Then, the cutoff distance dcut in some lo¬ 
cal region in (ap,ep,ip) is defined such that groups of 
N > A/min members cannot be produced with dcut by 
random fluctuations. All groups with N > A/min are 
then treated as meaningful asteroid families. Different re¬ 
searchers made different choices: A/min = 5 in Zappala 
et al. (1990), A/min = 100 in Parker et al. (2008), and 
A^min = 10-20 in most other publications. The two disad- 
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Fig. 2. — The Nysa-Polana complex, (a) The HCM applied to 
this region of the inner main belt reveals a major concentration of 
asteroids with 2.25 < ap < 2.48 AU and 0.13 < ep < 0.22. 
The shape of the concentration in the (ap,ep) projection is un¬ 
usual and difficult to interpret, (b) The WISE albedos of members 
of the Nysa-Polana complex: black for pp < 0.15 and gray for 
pv > 0.15. It becomes clear with the albedo information that 
the Nysa-Polana complex is two overlapping groups with distinct 
albedos. Furthermore, based on the V-shape criterion (Section 4), 
the low-albedo group is found to consist of two asteroid families 
(the Polana and Eulalia families; Walsh et al., 2013). The vertical 
feature at ap ~ 2.42 AU is the 1:2 mean motion resonance with 
Mars. 

vantages of this method are that: meaningful asteroid fam¬ 
ilies with N < A/min members are explicitly avoided, and 
<^cut(^min) depends on the population density in proper el¬ 
ement space and must be recomputed when a new classifi¬ 
cation is attempted from ever-growing catalogs. 

Another approach to this problem is to identify all 
groups, even if they have only a few members, and estab¬ 
lish their statistical significance a posteriori. Those that are 
judged to be insignificant are subsequently discarded and 
do not appear in the final lists. To determine the statistical 
significance of a group, one can generate mock distribu¬ 
tions and apply the HCM to them. For example, the high 
statistical significance of the Karin family, which is embed¬ 
ded in the much larger Koronis family, can be demonstrated 
by generating thousand orbital distributions corresponding 
to the Koronis family, and applying the HCM to each one 
(Nesvorny et al., 2002a). With dcut = 10ms“^,no con¬ 
centrations in this input can be found containing more than 
a few dozen members, while the Karin family currently has 
541 known members. Therefore, the Karin family is sig¬ 
nificant at a greater than the 99% level (again including a 
99% confidence interval of the estimate). A systematic ap¬ 
plication of this or similar statistical arguments can be quite 
laborious if many borderline cases need to be resolved. 

2.4 OVERLAP PROBLEM AND INTERLOPERS 

The overlap between different families has become more 
of a problem with a progressively larger share being taken 
in the proper element catalogs by small, km and sub-km 
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size asteroids. This is because small fragments are gener¬ 
ally launched at higher speeds, and are therefore initially 
spread in a larger volume in (ap,ep,ip) space. Mainly, 
however, the problem is caused by the larger mobility of 
small fragments due to the Yarkovsky effect. For example, 
the mean drift rate of a diameter D = 1 km main belt aster¬ 
oid is estimated to be 10“^ AU/m.y. (Bottke et al., 2006). 
The km-size members of a 1-g.y. old family are therefore 
expected to be dispersed over ^0.2 AU (the additional fac¬ 
tor of two accounts here for fragments having different spin 
orientations, and thus da/dt < 0 or da/dt > 0 ), which is 
roughly 1/5 of the extension of the whole main belt. In ad¬ 
dition, drifting asteroids encounter orbital resonances and 
can be dispersed by them in ep and ip as well. 

A good illustration of this is the case of the Flora and 
Vesta families in the inner main belt. To separate these 
families from each other down to their smallest members, 
the scope of the HCM can be restricted by an artificial cut 
in proper element space. Alternatively, one can first apply 
the HCM to the distribution of large members, thus identi¬ 
fying the core of each family, and then proceeding by trying 
to “attach” the small members to the core. This second step 
must use a lower dcut value than the first step to account for 
the denser population of smaller asteroids. In practice, this 
has been done by applying an absolute magnitude cutoff, 
H*, with H < Ff* for the core and H > iF* for the rest. 
In the low-i portion of the inner main belt, where the Flora 
and Vesta families reside, Milani et al. (2014) opted to use 
iF* = 15, and identified cores of families with Nmin = 17 
and dcut = 60 m s“^, and small members with Nmin = 42 
and dcut = 40 m s“^. 

Another solution to the overlap problem is to consider 
the physical properties of asteroids. Previously, the spec¬ 
troscopic observations of members of dynamical families 
have been used to: ( 1 ) establish the physical homogeneity 
of asteroid families (the difference between physical prop¬ 
erties of members of the same family tends to be smaller 
than the differences between physical properties of differ¬ 
ent families), and ( 2 ) identify large interlopers (asteroids 
classified as family members based on proper elements but 
having spectroscopic properties distinct from the bulk of the 
family). With the color and albedo data from the SDSS and 
WISE (Note 5), the physical homogeneity of asteroid fami¬ 
lies has been demonstrated to hold down to the smallest ob¬ 
servable members (Ivezic et al., 2001; Parker et al., 2008; 
see Figure [lb). A straightforward implication of this result 
is that the interior of each disrupted body was (relatively) 
homogeneous, at least on a scale comparable to the size of 
the observed fragments (^1-100 km) (Note 6 ). 

The physical homogeneity of asteroid families can be 
used to identify interlopers as those members of a dynam¬ 
ical family that have color and/or albedo significantly dis¬ 
tinct from the rest of the family. The number density of ap¬ 
parent color/albedo interlopers in a family can then be com¬ 
pared with the number density of the same color/albedo as¬ 
teroids in the immediate neighborhood of the family. Simi¬ 
lar densities are expected if the identified bodies are actual 


interlopers in the family. If, on the other hand, the density 
of color/albedo outliers in the family is found to be substan¬ 
tially higher than in the background, this may help to rule 
out the interloper premise, and instead indicate that: (i) the 
disrupted parent body may have been heterogeneous, or (ii) 
we are looking at two or more overlapping dynamical fam¬ 
ilies with distinct color/albedo properties. Finally, as for 
(ii), it is useful to verify whether the family members with 
different color/albedo properties also have different proper 
element distributions, as expected if breakups happened in 
two (slightly) different locations in proper element space 
(e.g., the Nysa-Polana complex; see Figure O. 

2.5 FAMILIES IN EXTENDED SPACE 

Another useful strategy is to include the color and/or 
albedo information directly in the clustering algorithm. 
This can be done by first separating the main belt into two 
(or more) populations according to their color and albedo 
properties. For example, asteroids in the S-complex can 
be separated from those in the C/X-complex based on the 
SDSS colors (Nesvorny et al. 2005), and the high-albedo 
asteroids can be separated from the low-albedo asteroids 
based on the albedo measurements of WISE (Masiero et 
al., 2013). The HCM is then applied to these popula¬ 
tions separately. This method is capable of identifying 
small/dispersed S-complex families in the C/X-type back¬ 
ground, and vice-versa, or low-albedo families in the high- 
albedo background, and vice-versa. It can also be useful to 
characterize the so-called family “halos” (Section 6.4). 

A more general method for including the color/albedo 
information in the clustering algorithm consists in the ap¬ 
plication of the HCM in space of increased dimension (e.g., 
Parker et al. 2008; Carruba et al., 2013a). When consider¬ 
ing the proper elements and SDSS colors, the distance in 5D 
can be defined as + k 2 {SC 2 )‘^), 

where d is the distance in 3D space of proper elements de¬ 
fined in Section 2.2, Ci and C 2 are two diagnostic colors 
defined from the SDSS (Ivezic et al., 2001; Nesvorny et 
al., 2005), and ki and F 2 are coefficients whose magni¬ 
tude is set to provide a good balance between the orbital 
and color dimensions (e.g., Nesvorny et al., 2006b). Sim¬ 
ilarly, we can define d\ = d? -\- n^a^kp{dpvY 4D) 
and ^4 = ^ 2 + in?a^kp{dpvY include the mea¬ 

surements of albedo pv from WISE. The 6/4 metric applies 
the strictest criteria on the family membership, because it 
requires that the family members have similar proper ele¬ 
ments, similar colors, and similar albedos. Note, however, 
that this metric can only be applied to a reduced set of main 
belt asteroids for which the proper elements, colors and 
albedos are simultaneously available (presently :^25,000; 
Figure dh). 

2.6 VERY YOUNG FAMILIES IN ORBITAL ELE¬ 
MENT SPACE 

Short after family’s creation, when the mutual gravity ef¬ 
fects among individual fragments cease to be important, the 
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Fig. 3.— Left panels: Proper elements of members of the Karin family. The size of each dark symbol is proportional to the diameter 
of a family member. Light gray dots indicate background bodies near the Koronis family. The black ellipses show the proper orbital 
elements of test bodies launched at 15 m s“^ from ap = 2.8661 AU, ep = 0.04449 and zp = 0.03692, assuming that / = 30° and 
a; + / = 45°, where / and u are the true anomaly and perihelion argument of the disrupted body at the time of the family-forming 
collision. Right panels; The convergence of angles at 5.8 m.y. ago demonstrates that the Karin family was created by a parent asteroid 
breakup at that time. The plot shows past orbital histories of ninety members of the Karin family: (top) the proper nodal longitude, and 
(bottom) the proper perihelion longitude. Values of these angles relative to (832) Karin are shown. At t = 5.8 m.y. (broken vertical 
line), the nodal longitudes and perihelion arguments of all ninety asteroids become nearly the same, as expected if these bodies had 
initially nearly the same orbits. Adapted from Nesvorny and Bottke (2004). 


fragments will separate from each other and find themselves 
moving on heliocentric orbits. Initially, they will have very 
tightly clustered orbits with nearly the same values of the 
osculating orbital angles O, w and A, where O is the nodal 
longitude, w is the apsidal longitude, and A is the mean lon¬ 
gitude. The debris cloud will be subsequently dispersed by 
the (i) Keplerian shear (different fragments are ejected with 
different velocity vectors, have slightly different values of 
the semimajor axis, and therefore different orbital periods) 
and (ii) differential precession of orbits produced by plane¬ 
tary perturbations. 

As for (i), the fragments will become fully dispersed 
along an orbit on a timescale = tt/ (adn/da) (Krb/ ^V) 
= {P/3){Vorh/^V), where P = 2-4 yr is the orbital period 
and SV is the ejection speed. With 6V = 1-100 m s“^, this 
gives Tn = 300-30,000 yr. Therefore, the dispersal of frag¬ 
ments along the orbit is relatively fast, and the clustering 
in A is not expected if a family is older than a few tens of 
thousand years. 

The dispersal of ft and w occurs on a time scale Tf = 
7r/{adf /da){Vorh/^V), where the frequency f = s or 
g. For example, ds/da = —70 arcsec yr“^ AU“^ and 
dgjda = 94 arcsec yr“^ AU“^ for the Karin family 
(a 2.865 AU). With 5V = 15 m s“^ (Nesvorny et al., 
2006a) and l^rb = 17.7 km s“^, this gives Tq = 3.8 m.y. 
and Tg = 2.8 m.y. Since fage > Ts and tage > Tg in this 


case, the distribution of ft and m for the Karin family is 
not expected to be clustered at the present time (Figure [3]) . 
Conversely, the clustering of ft and m would be expected 
for families with fage ^ 1 m.y. 

This expectation leads to the possibility that the families 
with fage ^ 1 m.y. could be detected in the catalogs of 
osculating orbital elements (Marsden, 1980; Bowell et al., 
1994), where they should show up as clusters in 5D space of 
a, e, i, and ft. The search in 5D space of the osculating 
orbital elements can be performed with the HCM method 
and metric d\ = d? ^ {na)‘^{kQ{6ft)‘^ + where 

d = (i(a, e, i) was defined in Section 2.2, and k^ and k^ 
new coefficients. [Different metric functions were studied 
by Rozek et al. (2011), who also pointed out that using the 
mean elements, instead of the osculating ones, can lead to 
more reliable results.] 

This method was first successfully used in practice for 
the identification of the Datura family (Nesvorny et al., 
2006c), and soon after for the discovery of the asteroid pairs 
(Vokrouhlicky and Nesvorny, 2008). The Datura family 
now consists of 15 known members ranging in size from 
:^10-km-diameter object (1270) Datura to sub-km frag¬ 
ments. They have ft and w clustered to within a few degrees 
near 98° and 357°, respectively. The age of the Datura fam¬ 
ily is only 530 ± 20 k.y., as estimated from the backward 
integration of orbits (Vokrouhlicky et al., 2009). Table 1 
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Semimajor axis ap (AU) 


Fig. 4. — The Veritas family. Here we plot the proper elements 
of 1294 members of the Veritas family identified from a new cat¬ 
alog (August 2014). (490) Veritas itself is marked by a trian¬ 
gle. The gray vertical strips indicate two three-body resonances 
which act to diffuse the orbits of family members in Cp. The 
black ellipse shows the orbital elements of test bodies launched 
at SV = 35 m s“^ from ap = 3.17 AU and ep = 0.0062, assum¬ 
ing that / = 30°, where / is the true anomaly of the disrupted 
body at the time of the family-forming collision. 

reports other notable cases of families with tage < 1 m.y. 

3. DETECTION OF RECENT BREAKUPS 

The detection of families with very young formation 
ages was one of the highlights of asteroid research in the 
past decade. A poster child of this exciting development is 
the Karin family, part of the larger Koronis family, that was 
shown to have formed only 5.8 ± 0.2 m.y. ago (Nesvorny 
et al., 2002a). The Karin family was identified by the tra¬ 
ditional means, using the HCM on proper elements. The 
diagonal shape of this family in the (ap, ep) projection is 
a telltale signature of a recent breakup, because this ini¬ 
tial shape is expected if a breakup occurs near the perihe¬ 
lion of the parent body orbit (Figure [S]). In this case, the 
tilt a of the family shape in the (ap, ep) projection is ex¬ 
pected to be, and indeed is in the case of the Karin family, 
tana = Aep/(Aap/ap) = (1 — ep) = 0.956. Thus, 
a 45°. All other known families, with an exception of 
the similarly-young Veritas family, have a ^ 0 (i.e., are 
nearly horizontal features in ap, ep). This is because these 
families are old and their original shape was stretched in ap 
by the Yarkovsky effect (e.g., DelTOro et al., 2004). 

The age of the Karin family has been established by nu¬ 
merically integrating the orbits of identified members back 
in time in an attempt to identify their past convergence (Fi¬ 
gure The past convergence is expected because the 
spread of the Karin family in ap, ep, ip indicates that the 
ejection speeds of observed fragments were only <15 m 
s“^. These low speeds imply a very tight initial distribu¬ 
tion (to within ^1°) of A, m and The backward inte¬ 
gration showed that the convergence occurred at 5.8 ± 0.2 


m.y. ago. In addition, the past convergence improved, with 
w and Q of all member orbits converging to within a de¬ 
gree, if the backward integration included the Yarkovsky 
drift (Nesvorny and Bottke, 2004). This was used to mea¬ 
sure the rate of the Yarkovsky drift for individual members 
of the family, determine their obliquities, and pin down the 
age of the Karin family to tage = 5.75 ± 0.05 m.y. 

The method of backward integration of orbits was ap¬ 
plied to several families (Nesvorny et al., 2002a, 2003, 
2008b; Novakovic et al., 2010, 2012a,b, 2014; Table [D 
some of these results will need to be verified). One of the 
interesting results that emerged from these studies is a pos¬ 
sible relationship between the young families and Main Belt 
Comets (MBCs; see chapter by Jewitt et al. in this volume). 
For example, (7968) P/133 Elst-Pizzaro can be linked to 
the Beagle family (tage ^ 10 m.y.; Nesvorny et al., 2008b), 
and (300163) P/2006 VW139 and P/2012 F5 (Gibbs) can 
be linked to small families that probably formed within the 
past 10 m.y. (Novakovic et al., 2012a, 2014). If this rela¬ 
tionship is confirmed by future studies, this can help us to 
understand how the MBCs are “activated”. 

Another notable case of a recent breakup is the Veritas 
family (Figure [H). It has previously been hypothesized that 
the Veritas family is <50 m.y. old (Milan! and Farinella 
1994). This claim was based on the argument that the 
largest member of the dynamical family, (490) Veritas, has 
chaotic dynamics and would be expected to diffuse away 
in ep from the rest of the family, if the family were older 
than ^50 m.y. [Note, however, that recent impact modeling 
may indicate that (490) Veritas is not a true member of the 
Veritas family (e.g., Michel et al., 2011).] A backward inte¬ 
gration of orbits confirmed the young age, and showed that 
the Veritas family formed only 8.3±0.5 m.y. ago (Nesvorny 
et al., 2003). 

A similarly young age was later obtained by an inde¬ 
pendent method, known as the “chaotic chronology”, based 
on tracking the evolution of orbits in one of the diffusive 
resonances that intersect the Veritas family (Tsiganis et al., 
2007, see also Knezevic and Pavlovic, 2002). Tsiganis et 
al. (2007) considered the chaotic diffusion of the Veri¬ 
tas family members in the 5J-2S-2 three-body resonance 
at 3.174 AU (Nesvorny and Morbidelli, 1999). Based on 
numerical integrations of chaotic orbits they estimated that 
the observed spread in the 5J-2S-2 resonance can be ob¬ 
tained for tage = 8.7 ± 1.7 m.y. Interestingly, the Veri¬ 
tas family is also intersected by the 7J-7S-2 resonance at 
3.168 AU. The observed distribution of eccentricities in 
this resonance is rather wide and cannot be explained by 
normal diffusion over the estimated age (the 7J-7S-2 res¬ 
onance is ^100 times less diffusive than the 5J-2S-2 res¬ 
onance). Perhaps the problem is with the HCM chaining, 
discussed in Section 2.2, which links unrelated asteroids in 
the 7J-7S-2 resonance, or the dynamical modeling is miss¬ 
ing some important ingredient. Novakovic et al. (2010) 
applied the method of chaotic chronology to the Theobalda 
family and found that the estimated age is consistent with 
that obtained from a backward integration of the regular or- 
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Table 1: Recently formed asteroid families. 


Eamily/Pair 

/age 

References 

Notes 

(832) Karin 

5.75 ± 0.05 m.y. 

Ne02,NB04 

2.1° dust band 

(158) Koronis(2) 

10-15 m.y. 

MH09 

near (832) Karin 

(490) Veritas 

8.3 ±0.5 m.y. 

Ne03,P06,T07 

9.3° band, late Miocene dust shower 

(656) Beagle 

^10 m.y. 

Ne08 

1.4° band, member Elst-Pizzaro? 

(778) Theobalda 

6.9 ± 2.3 m.y. 

NolO 

tage needs to be confirmed 

(1270) Datura 

530 ± 20 k.y. 

Ne06,V09 

identified in 5D, E/E dust band? 

(2384) Schulhof 

780 ± 100 k.y. 

VNll 

secondary breakup event? 

(4652) lannini 

^ 5 m.y. 

Ne03,W08 

chaotic dynamics 

(5438) Lorre 

1.9 ±0.3 m.y. 

Nol2a 

ip - 28° 

(14627) Emilkowalski 

220 ± 30 k.y. 

NV06 

only 3 members known 

(16598) 1992 YC2 

50-250 k.y. 

NV06 

only 3 members known 

(21509) Lucascavin 

300-800 k.y. 

NV06 

only 3 members known 

(300163) P/2006 VW139 

7.5 ±0.3 m.y. 

No 12b 

main belt comet 

P/2012 P5 (Gibbs) 

1.5 ±0.1 m.y. 

Nol4 

main belt comet 


References: NeOX = Nesvorny et al. (200X), NoXX = Novakovic et al. (20XX), VNXX = Vokrouhlicky and Nesvorny (20XX), VXX 
= Vokrouhlicky et al. (20XX), NB04 = Nesvorny and Bottke (2004), T07 = Tsiganis et al. (2007), F06 = Farley et al. (2006), W08 = 
Willman et al. (2008), NV06 = Nesvorny and Vokrouhlicky (2006), MH09 = Molnar and Haegert (2009). 


bits (tage = 6.9 ± 2.3 m.y.). 

There is a close relationship of the young asteroid fami¬ 
lies to the asteroid dust bands (see chapter by Jenniskens in 
this volume), which are strips of infrared emission running 
roughly parallel to the plane of the solar system (Low et al., 
1984). The three most prominent dust bands, known as a, 
/3 and 7, have previously been thought to originate in the 
Themis, Koronis and Eos families (Dermott et al., 1984). 
A detailed modeling, however, have shown that the sources 
of these dust bands are the recently-formed Karin, Beagle 
and Veritas families (Dermott et al., 2002; Nesvorny et al., 
2003, 2006, 2008b), mainly because: (1) the Veritas family 
with ip = 9.3° provides a better fit to the latitudinal posi¬ 
tion of the 7 band than the Eos family with ip 10°, and 
(2) the young families should now be more prolific sources 
of dust than the old families, because the dust production in 
a collisional cascade is expected to drop with time. A tracer 
of the Veritas family breakup has been found in measure¬ 
ments of extraterrestrial ^He in :^8.2-m.y.-old deep ocean 
sediments (Earley et al., 2006). 

4. V-SHAPE CRITERION 

A correct identification of the brightest/largest members 
in a family is important for several reasons. Eor example, 
the spectroscopic observations are magnitude limited and 
can typically only be conducted for bright targets. The in¬ 
terpretation of spectroscopic observations of asteroid fami¬ 
lies, and implications for the homogeneity/heterogeneity of 
their parent bodies, therefore depend on whether the bright 
asteroids are the actual members of a family, or not (e.g., 
Reddy et al., 2011). The large objects in families are also 
critically important for asteroid impact studies, in which 
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Eig. 5.— A recent application of the V-shape criterion to define 
the Polana and Eulalia families, and identify the largest members 
in these families (Walsh et al., 2013). The plot shows the absolute 
magnitude H of dark (pv < 0.1) asteroids with 0.1 < ep < 0.2 
and zp < 10° as a function of ap. The bottom panel illustrates 
the best-fit V-shape curves to the Polana (left) and Eulalia (right) 
families. (142) Polana and (495) Eulalia, the largest dark asteroids 
inside the left and right zones, respectively, are the likely largest 
members of the two families. 


the Size Erequency Distribution (SED) of family members 
is used to calibrate the results of impact experiments. Ei- 
nally, the uncertainty in the membership of the largest fam- 
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ily members critically affects our estimates of size of the 
disrupted parent bodies (e.g., Durda et al., 2007). 

While the interloper problem may not have an ideal solu¬ 
tion, there exists a straightforward method that can be used 
to remove obvious bright/large interlopers based on dynam¬ 
ical considerations. This method is inspired by the “V- 
shape” of asteroid families, which becomes apparent when 
the absolute magnitude H of family members is plotted 
against ap (Figure [5]). The V shape results from two pro¬ 
cesses. First, larger (smaller) fragments tend to be ejected 
at lower (higher) speeds, and thus tend to be located, on 
average, closer to (further away from) the family center 
(see chapter by Michel et al. in this volume). The sec¬ 
ond and typically more dominant ’V’ is contributed by the 
Yarkovsky effect. 

The Yarkovsky effect (YE) is a recoil force produced by 
anisotropic emission of thermal photons from an asteroid 
surface (see chapter by Vokrouhlicky et al. in this volume). 
The diurnal component of the YE, which is more impor¬ 
tant for asteroid-sized bodies than its seasonal counterpart, 
can increase the semimajor axis of an asteroid with a pro¬ 
grade rotation and decrease the semimajor axis for a retro¬ 
grade rotation. The semimajor axis drift is generally given 
by dap/dt = const, cos 0/{apD), where 0 is the asteroid 
obliquity, D is the effective diameter, and the constant de¬ 
pends on material properties. The maximum drift occurs 
for 0 = 0° or 0 = 180°. Thus, the envelope of the dis¬ 
tribution of family members in (ap, 77) is expected to fol¬ 
low |ap — ad = CyeIO^/^, where ac is the family cen¬ 
ter (often assumed to coincide with the largest fragment), 
and Cye is a constant related to tage (Section 5). Interest¬ 
ingly, if the ejection speed SV oc 1/77, as found for the 
youngest families (e.g., Nesvorny et al., 2006a) and some 
laboratory experiments (e.g., Fujiwara et al., 1989), then the 
ejection velocity field will produce the same dependence, 
|<^p — del = CevIO^/^, where Cev is related to the mag¬ 
nitude of the ejection velocities. 

These considerations allow us to define the V-shape cri¬ 
terion (e.g., Nesvorny et al., 2003). Consider a family ex¬ 
tracted by the HCM as described in Section 2.2. Fit an enve¬ 
lope to the distribution of small family members in (ap, 77) 
using the functional dependence between ap and 77 defined 
above. The envelope is then continued to low 77 values, and 
the bright members of the HCM family that fall outside the 
envelope boundaries can be rejected. On the other hand, 
the brightest HCM family members that fall within the en¬ 
velope boundaries are good candidates for the largest family 
members. This method is illustrated in Figure [3 

In practice, a number of additional effects can com¬ 
plicate the application of the V-shape criterion described 
above. For example, strong nearby resonances can remove 
family members that drifted into them, thus producing cut¬ 
offs of the ap — ac distribution beyond which no family 
members can be identified (e.g., the 3:1 resonance in Fig¬ 
ure B- Also, family members can be displaced in ap by 
encounters with (1) Ceres and other massive main belt as¬ 
teroids (Nesvorny et al., 2002b; Carruba et al., 2003,2007a, 


2012, 2013b; Delisle and Laskar, 2012). In addition, the 
physics of large scale collisions is still poorly understood, 
and the possibility that some large fragments can be accel¬ 
erated to very high speeds cannot be ruled out. Therefore, 
while the V-shape criterion defined above is a useful guide, 
it cannot be rigidly applied. 

Possibly the best way to deal with this issue is to de¬ 
fine the value of Co = 10“^/^|ap — ad that best fits the 
V-shaped family envelope and report Cj/Co, where Cj = 
10“^i/^|ap^j — ad, for each family member j identified 
by the HCM. Asteroids with |Cj/Co| > 1 can be flagged 
(but not removed), because they are potential interlopers in 
a dynamical family according to the V-shape criterion. The 
results can then be cross-linked with the spectroscopic data 
(or colors, or albedo measurements) to determine whether 
there is a good correspondence between the flagged bodies 
and (suspected) spectroscopic interlopers. An illustration of 
this procedure can be found in Vokrouhlicky et al. (2006a) 
who examined the Eos family. They found that many large 
bodies in the HCM family with | Cj /Co | > 1 have physi¬ 
cal properties that are indeed incompatible (mainly dark C 
types) with the bulk of the Eos family (mainly brighter K 
types). 

5. FAMILY AGE ESTIMATION 

The method of backward integrations of orbits described 
in Section 3 cannot be used to determine ages of the fam¬ 
ilies much older than ^10 m.y., mainly because the or¬ 
bital evolution of main belt asteroids is generally unpre¬ 
dictable on long timescales and sensitively depends on non- 
gravitational effects that are difficult to model with the 
needed precision. Instead, the age of an old family can 
be estimated by a statistical method, which is based on the 
general notion that the spread of an asteroid family in ap 
increases over time as its members drift away due to the 
Yarkovsky effect. Expressed in terms of the equations dis¬ 
cussed in the previous section, the age can be estimated as 



where p is the asteroid bulk density and pv is the visual 
albedo. While the equation above is scaled to typical val¬ 
ues expected for an S-type asteroid, a change to p = 1.5 g 
cm“^ and pv = 0.05, which would be more appropriate for 
a dark C-type asteroid, produces two multiplication factors 
that nearly compensate each other. The biggest uncertainty 
in the inversion from Co to tage lies in the unknown density 
factor. Additional uncertainty, not explicitly apparent from 
Eq. (2), stems from the dependence of the drift rate on sur¬ 
face conductivity K. Together, the conversion from Co to 
tage has an uncertainty of about a factor of 2. 

More fundamentally, Eq. ([T]) neglects complicating fac¬ 
tors such as the contribution of the original ejection field 
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and nearby resonances that can remove drifting bodies. If 
C = Cye + C'ev, under the assumptions on the ejection 
speeds discussed in Section 4, the two effects cannot be de¬ 
coupled from each other, and tage estimated from Cq will 
always overshoot the real age of the family. This happens, 
in essence, because this simple method only fits the family 
envelope. 

Vokrouhlicky et al. (2006a,b) have developed a more 
general statistical method that uses the actual distribution 
of family members within the V-shape envelope. It works 
as follows. First, the code generates a new-born family. The 
distribution of the ejection speeds of fragments is approx¬ 
imated by the Gaussian distribution in each velocity com¬ 
ponent with the size-dependent standard deviation 5V = 
1 / 5(5 km/D), where V 5 is a free parameter. This functional 
dependence generally provides a good match to the ejection 
speeds inferred from the young families and impact simu¬ 
lations (Michel et al., 2001, 2003, 2004; Nesvorny et al., 
2002a, 2006a; Durda et al., 2004, 2007). The Yarkovsky- 
YORP (hereafter YY) code then evolves the orbits of frag¬ 
ments accounting for the semimajor axis drift due to the 
YE, and the spin evolution from the YORP effect (Rubin- 
cam, 2000). To speed up the calculation, planetary pertur¬ 
bations were not taken into account in the original YY code 
(but see Masiero et al., 2012; results based on full A^-body 
integrations are discussed in the following section). 

When applied to an asteroid family, the goal is to find a 
combination of parameters V 5 and tage that best fits the ob¬ 
served (ap, 77) distribution in the family. Since the model 
does not contain an SFD-related evolution component, this 
can be conveniently done by simply fitting the observed dis¬ 
tribution of Cj . The results are found to be credible in cases 
such as the Erigone family (Figure©, where the YY model 
is capable of adequately representing the observed C distri¬ 
bution, which shows void space near C = 0, a maximum at 
intermediate values of C, and a relatively sharp drop toward 
the largest C values seen in the family. This dependence is 
produced by the YORP effect which tilts the spin axis away 
from the orbital plane and therefore, through the cos 0 de¬ 
pendence of the YE, tends to maximize the drift rates. This 
gives some asteroid families the characteristic appearance 
in (ap, 77), which different authors, with noticeably differ¬ 
ent gifts for subtleties of poetic expression, called “ears”, 
“wings” or “petals”. 

But not all families have ears, wings or petals. In fact, 
most families have more uniform distribution of C that goes 
all the way from 0 to Cmax- This is thought to have some¬ 
thing to do with how the YE and YORP operate on g.y.- 
long timescales. It is perhaps related to the variability of 
the spin-related YORP torque, which sensitively depends 
on small perturbations of the asteroid surface. As a result 
of this dependence, the spin rate of each individual aster¬ 
oid can undergo a random walk (stochastic YORP; Bottke 
et al., 2015) and, when spun down completely, the aster¬ 
oid spin axis would chaotically tumble for some time, and 
subsequently reorient. This process could mix up the total 
drifts suffered by family members and, rather unpoetically. 


remove the ears, wings or petals from the old families. 

The application of the standard YY code (and, for that 
matter, also of the simple method based on the ‘V’-shape 
envelope) is problematic in cases when a family is cut by 
nearby mean motion resonances and looks like a box in 
(ap, 77), or when Cev is expected to be larger than Cye, 
perhaps because the family’s parent was an object with 
large escape speed (e.g., the Vesta family; (4) Vesta has 
Vesc — 360 m s“^), or because the family formed recently 
(e.g., the Datura, Karin, Veritas families have Cye — 0). 
A complete list of families to which the YY code was suc¬ 
cessfully applied so far is: 

• Agnia (fage = 100 ± 100 m.y., V 5 15 m s“^) 

• Massalia (fage = 150 ± 50 m.y., V 5 20 m s“^) 

• Baptistina (tage = 160 ± 50 m.y., V 5 40 m s“^) 

• Merxia (fage = 250 ± 100 m.y., V 5 25 m s“^) 

• Astrid (fage = 250 ± 100 m.y., V 5 15 m s“^) 

• Erigone (tage = 300 ± 100 m.y., V 5 30 m s“^) 

• Eos (fage = 1-3 ± 0.5 g.y., V 5 70 m s“^) 

• Tina (fage = 170 ± 50 g.y., V 5 20 m s“^). 

Here we have taken the liberty to update and round off 
the estimates, and give more generous errors than in the 
original publications (Vokrouhlicky et al., 2006a,b,c; Bot¬ 
tke et al., 2007; Carruba and Morbidelli, 2011; see also Car- 
ruba, 2009a for the Padua family). Note that these errors do 
not include the uncertainty of about a factor of 2 from the 
poorly known bulk density and surface conductivity of the 
asteroids in question. Including this uncertainty, Masiero et 
al. (2012) found that the best-fitting age of the Baptistina 
family can be anywhere between 140 and 320 m.y. 

The estimated ejection speeds are V 5 = 15-50 m s“^, 
except for the Eos family, which formed in a breakup of 
a very large parent asteroid (77 pb ^ 300 km). These re¬ 
sults are consistent with the ejection speeds inferred from 
the young Karin family, which has V 5 15ms“^ for a 
relatively small parent body (79 pb — 35 km; Nesvorny et 
al., 2006a). The ejection speeds contribute by :^20% (for 
oldest Eos) to 50% (for youngest Agnia) to the total family 
spread in the semimajor axis. Ignoring this contribution, as 
in Eq. (2), would thus lead to an overestimate of tage by 
^20-50%. While one must therefore be careful in applying 
Eq. (2) to the small/young families that did not have enough 
time to significantly spread by the YE, the effect of the ejec¬ 
tion speeds should be less of an issue for old families. 

6. DYNAMICAL EVOLUTION 

6.1 INITIAL STATE 

The dynamical evolution of asteroids in families is sim¬ 
ilar to the dynamical evolution of main belt asteroids in 
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Fig. 6 .— (a) The Erigone family projected onto a plane of proper semimajor axis versus absolute magnitude H\ (163) Erigone is the 
filled star. The family has been separated into two clouds (up < 2.37 AU; Up > 2.37 AU) by the Yarkovsky/YORP evolution, (b) A 
comparison between model results (solid line) and binned Erigone family (gray dots; see Vokrouhlicky et al., 2006b). The error bars are 
the square root of the number of bodies in each bin. The a;-axis is the distance of family members from the the family center. Based on 
this result, Vokrouhlicky et al. (2006b) estimated that tage = 28 OI 5 Q m s“^, where the error bars do not include 

the uncertainty originating from uncertain material properties (e.g., density, surface conductivity). 


general. Studying the dynamical evolution of individual 
families is useful in this context, because we more or less 
know how the families should look like initially. Things 
may thus be learned by comparing these ideal initial states 
with how different families look now, after having dynami¬ 
cally evolved since their formation. The dynamical studies 
can also often provide an independent estimate of tage- 

Assuming that SV C V^rb. the initial shape of fam¬ 
ilies in can be obtained from the Gauss equa¬ 

tions (e.g., Zappala et al., 2002), which map the initial 
velocity perturbation = (Vr^Vt^Vz), where Ur, Vt 
and Vz are the radial, tangential and vertical components 
of the velocity vector, to the change of orbital elements 
SFi = {6a, 6e, 6i). If the ejection velocity field is (roughly) 
isotropic, the Gauss equations imply that initial families 
should (roughly) be ellipsoids in (a, e, i) centered at the ref¬ 
erence orbit (a*, e*, i*). The transformation from (a, e, i) 
to (ap,ep,ip) preserves the shape, but maps 
onto (up, Cp, ip) such that, in general, Up 7 ^ a*, Cp 7 ^ e* 
andi^ 7^ i*. 

The shape of the ellipsoids in (ap, ep, ip) is controlled 
by the true anomaly / and the argument of perihelion uj of 
the parent body at the time of the family-forming breakup. 
The projected distribution onto the (ap, ep) plane is a tilted 
ellipse with tightly correlated ap and ep if the breakup hap¬ 
pened near perihelion (see Figure [3^) or tightly anticorre¬ 
lated ap and ep if the breakup happened near aphelion. The 
two recently-formed families for which this shape is clearly 
discernible, the Karin and Veritas families, have correlated 
ap and ep, implying that |/| — 30° (Figures [5] and ID). 

The projected initial distribution onto the (ap, ip) plane 
is an ellipse with horizontal long axis and vertical short 
axis. The short-to-long axis ratio is roughly given by 
cos(cc; + /)Vz/Vt. Thus, breakups near the ascending 


{uj-\- f c^O) and descending {uj-\- f ir) nodes should pro¬ 
duce ‘fat’ ellipses while those with cj + / = zb 7 r /2 should 
make ‘squashed’ ellipses with Sip 0. While the Karin 
family neatly fits in this framework with uj-\-f it (Figu¬ 
re O, the Veritas family shows large Sip values, indicating 
that the ejection velocity field should have been anisotropic 
with Vz some 2-4 times larger than Ur- 

The reference orbit (ap, Cp, ip) is often taken to coin¬ 
cide with the proper orbit of a largest family member. This 
should be fine for families produced in cratering or mildly 
catastrophic events, where the orbital elements of the im¬ 
pacted body presumably did not change much by the im¬ 
pact. For the catastrophic and highly catastrophic breakups, 
however, the largest surviving remnant is relatively small 
and can be significantly displaced from the family’s cen¬ 
ter. For example, (832) Karin, the largest :^17-km-diameter 
member of the Karin family produced by a catastrophic 
breakup of a ^40-km-diameter parent body (mass ratio 
^0.08; Nesvorny et al., 2006a), is displaced by —0.002 AU 
from the family center (:^ 20 % of the whole extension of the 
Karin family in ap). This shows that, in general, the posi¬ 
tion of the largest fragment does not need to perfectly coin¬ 
cide with the family center, and has implication for the V- 
shape criterion discussed in Section 4 (where an allowance 
needs to be given for a possible displacement). 

6.2 DYNAMICS ON GIGAYEAR TIMESCALES 

An overwhelming majority of the observed asteroid fam¬ 
ilies are not simple Gaussian ellipsoids. While this was not 
fully appreciated at the time of the Asteroids III book, to¬ 
day’s perspective on this issue is clear: the families were 
stretched in ap as their members drifted away from their 
original orbits by the Yarkovsky effect. The asteroid fami¬ 
lies found in the present main belt are therefore nearly hor- 
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Fig. 7.— Dynamical evolution of the Eos family. From left to right, the panels show the: (a) observed family and its halo in the (ap, ep) 
projection, (b) assumed initial shape of the family, and (c) family’s structure after 1.7 g.y. In (a), we plot all asteroids with Eos-family 
colors (0.0 < a* <0.1 mag and —0.03 < i — z < 0.08 mag; see Ivezic et al. (2001) for the definition of color indices from the SDSS). 
The size of a symbol is inversely proportional to absolute magnitude H. The boxes approximately delimit the extent of the core and 
halo of the Eos family. In (b), 6545 test particles were distributed with assumed isotropic ejection velocities, V 5 = 93ms“^,/ = 150° 
and (jj = 30°. Nearly all initial particles fall within the family core. In (c), an A^-body integrator was used to dynamically evolve the 
orbits of the test particles over 1.7 g.y. The integration included gravitational perturbations from planets, and the Yarkovsky and YORP 
effects. The vertical lines show the locations of several resonances that contributed to spreading of orbits in ep (7:3, 9:4 and 11:5 with 
Jupiter, also 3J-2S-1 and zi = g s — ge — sq = 0). Adapted from Broz and Morbidelli (2013). 


izontal and elongated structures in (ap,ep) and (ap,ip). 
This shows that the original ejection velocity field can¬ 
not be easily reconstructed by simply mapping back to¬ 
day’s (ap, ep, ip) to (1 /r, Vt^Vz) from the Gauss equations 
(Note 7). 

Moreover, many asteroid families have weird shapes 
which, taken at the face value, would imply funny and 
clearly implausible ejection velocity fields. A prime exam¬ 
ple of this, briefly mentioned in section 2.2, is the Koronis 
family (Bottke et al., 2001). Since the case of the Koronis 
family was covered in the Asteroids III book (chapter by 
Bottke et al., 2002), we do not discuss it here. Instead, we 
concentrate on the results of new dynamical studies, many 
of which have been inspired by the Koronis family case. 
The dynamical effects found in these studies fall into three 
broad categories: 

(1) Members drifting in ap encounter a mean motion 
resonance with one of the planets (mainly Jupiter, Mars or 
Earth; see chapter Nesvorny et al., 2002c in Asteroids III 
book). If the resonance is strong enough (e.g., 3:1, 2:1, 
or 5:2 with Jupiter), the orbit will chaotically wander near 
the resonance border, its eccentricity will subsequently in¬ 
crease, and the body will be removed from the main belt and 
transferred onto a planet-crossing orbit (Wisdom, 1982). If 
the resonance is weak, or if the asteroid is small and drifts 
fast in ap, the orbit can cross the resonance, perhaps suf¬ 
fering a discontinuity in ep during the crossing, and will 
continue drifting on the other side. If the resonance is weak 
and the drift rate is not too large, the orbit can be captured 
in the resonance and will slowly diffuse to larger or smaller 
eccentricities. It may later be released from the resonance 
with ep that can be substantially different from the origi¬ 


nal value. The effects of mean motion resonances on ip are 
generally smaller, because the eccentricity terms tend to be 
more important in the resonant potential. The inclination 
terms are important for orbits with ip > 10°. A good ex¬ 
ample of this is the Pallas family with ip 33° (Carruba et 
al., 2011 ). 

(2) Drifting members meet a secular resonance. The 
secular resonances are located along curved manifolds in 
(ap,ep,ip) space (Knezevic et al., 1991). Depending on 
the type and local curvature of the secular resonance, and 
asteroid’s da/dt, the orbit can be trapped inside the reso¬ 
nance and start sliding along it, or it can cross the resonance 
with a noticeably large change of ep and/or ip. A good ex¬ 
ample of the former case are orbits in the Eos family sliding 
along the ^1 = ^ — ^6 + “ <^6 = 0 resonance (Vokrouhlicky 
et al., 2006a). An example of the latter case is the Koronis 
family, where eccentricities change as a result of crossing of 
the g -b 2^5 — 3^6 = 0 resonance (Bottke et al., 2001). If the 
secular resonance in question only includes the g (or s) fre¬ 
quency, effects on ep (or ip) are expected. If the resonance 
includes both the g and s frequencies, both ep and ip can 
be effected. If the orbit is captured in a resonance with the 
g and s frequencies, it will slide along the local gradient of 
the resonant manifold with changes of ep and ip depending 
on the local geometry. 

(3) Encounters with (1) Ceres and other massive aster¬ 
oids produce additional changes of ap, ep and ip (Nesvorny 
et al., 2002b; Carruba et al., 2003, 2007a, 2012, 2013b; 
Delisle and Laskar, 2012). These changes are typically 
smaller than those from the Yarkovsky effect on ap and res¬ 
onances on ep and ip. They are, however, not negligible. 
The effect of encounters can be approximated by a random 
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walk (see Carmba et al., 2007a for a discussion). The mean 
changes of ap, ep and ip increase with time roughly as ^/i. 
The asteroid families become puffed out as a result of en¬ 
counters, and faster so initially than at later times, because 
of the nature of the random walk. Also, a small fraction 
of family members, in some cases perhaps including the 
largest remnant, can have their orbits substantially affected 
by a rare, very close encounter. Additional perturbations of 
asteroid orbits arise from the linear momentum transfer dur¬ 
ing non-disruptive collisions (Dell’Oro and Cellino, 2007). 

6.3 DISCUSSION OF DYNAMICAL STUDIES 

The Koronis family case (Bottke et al., 2001) sparked 
much interest in studies of the dynamical evolution of as¬ 
teroid families on very long timescales. Here we review 
several of these studies roughly in the chronological order. 
The goal of this text is to illustrate the dynamical processes 
discussed in the previous section on specific cases. 

Nesvorny et al. (2002b) considered the dynamical evolu¬ 
tion of the Flora family. The Flora family is located near the 
inner border of the main belt, where numerous mean motion 
resonances with Mars and Earth produce slow diffusion of 
Cp and ip. The numerical integration of orbits showed how 
the overall extent of the Flora family in ep and ip increases 
with time. The present width of the Flora family in ep and 
ip was obtained in this study after t 0.5 g.y. even if 
the initial distribution of fragments in ep and ip was very 
tight. The Flora family expansion saturates for t > 0.5 g.y., 
because the Flora family members that diffused to large ec¬ 
centricities are removed from the main belt by encounters 
with Mars (the Flora family is an important source of chon- 
dritic near-Earth asteroids (NBAs); Vernazza et al., 2008). 
The present spread of the Flora family in ap, mainly con¬ 
tributed by the Yarkovsky effect, indicates tage ^ 1 g-Y- 

Carruba et al. (2005) studied the dynamical evolution 
of the Vesta family. The main motivation for this study 
was the fact that several inner main belt asteroids, such as 
(956) Elisa and (809) Lundia, have been classified as V- 
types from previous spectroscopic observations (Florczak et 
al., 2002), indicating that they may be pieces of the basaltic 
crust of (4) Vesta. These asteroids, however, have orbits 
rather distant form that of (4) Vesta and are not members of 
the Vesta’s dynamical family even if a very large cutoff dis¬ 
tance is used. It was therefore presumed that they: (i) have 
dynamically evolved to their current orbits from the Vesta 
family, or (ii) are pieces of differentiated asteroids unrelated 
to (4) Vesta. Carruba et al. (2005) found that the interplay 
of the Yarkovsky drift and the Z 2 = 2{g — qq) ^ s — sq = 0 
resonance produces complex dynamical behavior that can, 
indeed, explain the orbits of (956) Elisa and (809) Lundia, 
assuming that the Vesta family is at least ^1 g.y. old. This 
gives support to (i). 

In a follow-up study, Nesvorny et al. (2008a) performed 
a numerical integration of 6,600 Vesta fragments over 2 g.y. 
They found that most V-type asteroids in the inner main belt 
can be explained by being ejected from (4) Vesta and dy¬ 


namically evolving to their current orbits outside the Vesta 
family. These V-type “fugitives” have been used to con¬ 
strain the age of the Vesta family, consistently with find¬ 
ings of Carruba et al. (2005), to tage ^ 1 g-Y- Since pre¬ 
vious collisional modeling of the Vesta family suggested 
tage ^ 1 g-Y- (Marzari et al., 1999), the most likely age of 
the Vesta family that can be inferred from these studies is 
^age ^ 1 g-Y- This agrees well with the age of the :^500-km- 
diameter Rheasilvia basin on (4) Vesta inferred from crater 
counts (:^ 1 g.y.; Marchi et al., 2012) (Note 8). 

Vokrouhlicky et al. (2006a) studied the dynamical evo¬ 
lution of the Eos family. The Eos family has a compli¬ 
cated structure in proper element space leading some au¬ 
thors to divide it in several distinct families (e.g., Milani 
et al., 2014). Diagnostically, however, the Eos family, al¬ 
though somewhat physically heterogeneous, has the color, 
albedo and spectral properties that contrast with the local, 
predominantly C-type background in the outer asteroid belt. 
This suggests that this is a single family. As we discuss be¬ 
low, the complicated structure of the Eos family arises from 
of the presence of several mean motion and secular reso¬ 
nances. 

To start with, Vokrouhlicky et al. (2006a) showed that 
the Eos family members drifting by the Yarkovsky effect 
into the 7:3 resonance with Jupiter are removed (see Eig- 
ure[7]). This cuts the family at 2.957 AU. Members drifting 
with da/dt > 0, on the other hand, will encounter the 9:4 
resonance at 3.03 AU. This resonance, being of higher order 
and thus weaker, is not an unpenetrable barrier, especially 
for smaller members with higher drift rates. The estimated 
fraction of bodies that can cross the 9:4 resonance is <10% 
for 77 < 12 but reaches :^35%for77 = 16. This is con¬ 
sistent with the magnitude distributions of the Eos family 
members on both sides of the 9:4 resonance. The larger 
dispersion of the part of the Eos family with ap > 3.03 AU 
is contributed by perturbations of ep and ip during the 9:4 
resonance crossing (Eigure[7]). Einally, many orbits in the 
central part of the Eos family are trapped in the secular res¬ 
onance zi = g ^ s — ge — sq = 0, and slide along it while 
drifting in ap (Note 9). 

Einally, we discuss additional processes whose signifi¬ 
cance is shadowed by the Yarkovsky effect and resonances, 
but which can be important in some cases. Nesvorny et 
al. (2002b) considered encounters with (1) Ceres and found 
that the characteristic change of the semimajor axis due to 
these encounters is Aa 0.001 AU over 100 m.y. Assum¬ 
ing that the scattering effect of encounters can be described 
by a random walk with Aa oc \/t, the expected changes 
over 1 g.y. and 4 g.y. are :^0.003 AU and :^0.007 AU, re¬ 
spectively. The orbital changes were found to be larger for 
orbits similar to that of (1) Ceres, because the orbital prox¬ 
imity leads to lower encounter speeds, and larger gravita¬ 
tional perturbations during the low-speed encounters. Car¬ 
ruba et al. (2003) studied the effect of encounters on the 
Adeona and Gefion families, both located near (1) Ceres in 
proper element space. They found that the semimajor axis 
of members of the Adeona and Gefion families can change 
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by up to ^0.01 AU over the estimated age of these families. 

With similar motivation, Carruba et al. (2007a) con¬ 
sidered the effect of encounters of the Vesta family mem¬ 
bers with (4) Vesta. They found the characteristic changes 
Aa = 0.002 AU, Ae = 0.002 and Ai = 0.06° over 
100 m.y. In a follow-up work. Delisle and Laskar (2012) in¬ 
cluded the effects of eleven largest asteroids. They showed 
that encounters of the Vesta family members with (4) Vesta 
and (1) Ceres are dominant, contributing roughly by 64% 
and 36% to the total changes, respectively. The functional 
dependence Aa = 1.6 x 10“^^t/l m.y. AU was used in 
this work to extrapolate the results to longer time intervals. 
Moreover, Carruba et al. (2013b) studied the influence of 
these effects on the Pallas, Hygiea and Euphrosyne fami¬ 
lies. They showed that the effects of (2) Pallas -the 3rd 
most massive main-belt asteroid- on the Pallas family are 
very small, because these asteroids have high orbital in¬ 
clinations (ip 33°), lower frequency of encounters and 
higher-than-average encounter speeds. 

DeirOro and Cellino (2007) pointed out that orbits of 
main belt asteroids can change as a result of the linear mo¬ 
mentum transfer during non-destructive collisions. They 
found that the expected semimajor axis change from these 
collisions for a D = 50-km main-belt asteroid is Aa ^ 
10“^ AU over 100 m.y. (with the scaling laws from Benz 
and Asphaug, 1999). This is an order of magnitude lower 
than the change expected from close encounters with large 
asteroids and comparable to the sluggish drift rate expected 
from the Yarkovsky effect for D = 50 km. For D < 50 km, 
the orbital changes from non-destructive collisions sensi¬ 
tively depend on several unknown parameters, such as the 
SFD of sub-km main-belt asteroids, but the general trend 
is such that Aa drops with decreasing D (assuming that 
the cumulative SFD index is <4; DelTOro and Cellino, 
2007). Since Aa is independent of D for encounters with 
(1) Ceres, and Aa oc 1/D for the Yarkovsky force, these 
two effects outrun the non-destructive collisions for D < 
50 km. This limits the significance of non-destructive colli¬ 
sions for the dynamical evolution of asteroid families. Their 
effect on ep and ip is also minor. 

6.4 FAMILY HALOS 

When dynamical families are identified and removed 
from the main belt, they leave behind holes in the distri¬ 
bution of proper elements that are surrounded by regions 
with increased asteroid density (Figure [8]). These peripheral 
regions are known as the family halos. The families and 
family halos surrounding them are clearly related, which 
can most conveniently be demonstrated by considering their 
physical properties. For example, the Koronis family and its 
halo consist of bright asteroids (mean albedo pv = 0.15) 
that are classified as S in the asteroid taxonomy (moder¬ 
ate spectral slope and shallow absorption band near 1 pm). 
These properties contrast with the local background in the 
outer main belt, which is mostly dark (pv — 0.05) and C- 
type (featureless neutral spectrum). The Eos family and its 



Fig. 8.— The halos of families in the inner main belt (2 < ap < 
2.5 AU). The dynamical families were identified by the HCM in 
3D of proper elements and were removed, leaving behind what 
resembles a fuzzy chalk contour of victim’s body removed from 
a crime scene. The head, body and bottom leg correspond to the 
locations of the Vesta, Flora and Nysa-Polana families. The halos 
surrounding these families, mainly composed of small asteroids, 
are rather dramatic in this projection. This figure illustrates a ma¬ 
jor weakness of the traditional family-identification method based 
on the HCM in 3D space of proper elements. 


halo (K-type, pv = 0.13) also stand out from the dark outer 
belt background. 

The distinction between families and family halos is, in 
a sense, formal, because it appears as a consequence of the 
identification method. Indeed, to extract the family halo 
by applying the HCM in 3D space of proper elements, one 
would need to increase the value of the d^ut parameter. If 
it is increased, however, the identification would fail due 
to the problem of chaining (Section 2.2). This is a conse¬ 
quence of the fact that the number density of halo asteroids 
is comparable to, or only slightly larger than, the number 
density of background asteroids. A better method for ex¬ 
tracting the families jointly with their halos consists in ap¬ 
plying the HCM in space of extended dimensions, where 
physical properties are used in addition to the proper el¬ 
ements (Section 2.5; Parker et al., 2008; Carruba et al., 
2013a; Masiero et al., 2013). 

The origin of family halos is contributed by at least two 
processes. In the case of the Eos family, the observed dis¬ 
persion of halo in ep and ip is too large (0-0.15 and 7- 
14°, respectively) to be related to the the ejection speeds. 
Instead, Broz and Morbidelli (2013) showed how the Eos 
family halo gradually appears as a result of complex orbital 
dynamics in the Eos family region (Figure [7]). The age of 
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the Eos family inferred from that work, 1.5 g.y. < tage < 
2.2 g.y., is larger than the one obtained from the YY code 
(tage = 1-3 ± 0.2 g.y.; Vokrouhlicky et al., 2006a and Sec¬ 
tion 5), but the two estimates are consistent at the la level. 

The Koronis family and its halo, on the other hand, 
are located in a dynamically quiet region of the main belt 
(2.82 < ap < 2.95 AU, ep 0.045, ip 2°) where 
no major resonances exist (except for ^ + 2^5 — Sge = 0 ; 
Bottke et al., 2001). The long-term dynamics of the Koro¬ 
nis family therefore has a relatively small influence on the 
overall width of the family in ep and ip. In contrast, the 
present shift of the Koronis family/halo members from the 
family’s center is relatively large (up to 0.015 in ep and 
0.5° in ip). This can plausibly be explained by the ejection 
speeds. 

Assuming an isotropic ejection velocity held and 6V = 
Vs{b km/ D) with V5 = 50 m s“^, we And that 6ap/ap 
25V/Vorh — 0.018 for T) = 5 km fragments. This repre¬ 
sents a small share of the Koronis family width in ap (the 
population of D = b km members stretches from the 5:2 
resonance at 2.82 AU to the 7:3 resonance at 2.95 AU), 
consistently with the idea that the observed spread in ap 
is mainly due to the Yarkovsky effect. With D = 1.5 km, 
which is roughly the size of the smallest asteroids observed 
in this part of the main belt, the expected changes are ^ep = 
(1-2) X 5V/Vorh ^ 0.015 and 5ip < 5V/Vorh ^ 0.6°. 
This is comparable to the present width of the Koronis fam¬ 
ily/halo in ep and ip 

Another notable argument in this context is that the Ko¬ 
ronis halo is populated by small asteroids with diameters 
near the minimum detectable size (while the big ones are 
inside the family core). This suggest a certain size de¬ 
pendency of the process that created the halo, presumably 
related to the size dependency of the ejection speeds (see 
chapter by Michel et al. in this volume). This important 
trend was pointed out by Cellino et al. (2004), and was 
further discussed in Cellino et al. (2009). 

7. SYNTHESIS 

7.1 METHOD, FAMILY NAMES AND FIN 

Here we present a synthesis of asteroid families ex¬ 
tracted from many past publications (Mothe-Diniz et al., 
2005; Nesvorny et al., 2005; Gil-Hutton, 2006; Parker et 
al., 2008; Nesvorny 2010, 2012; Novakovic et al., 2011; 
Broz et al., 2013; Masiero et al., 2011, 2013; Carruba et 
al., 2013a; Milani et al., 2014). The synthesis is based on 
the numerically computed proper elements (AstDyS cata¬ 
log from August 2014 with the proper elements for 384,337 
numbered asteroids; the proper elements of 4016 Jupiter 
Trojans were taken from Rozehnal and Broz 2013) and the 
latest releases of data from the SDSS and WISE. 

Eirst, we collected a long list of all (they were 228) fam¬ 
ilies reported in these publications. Second, we checked to 
see whether the clustering in proper elements of each re¬ 
ported group is meaningful, and whether the group has rea¬ 


sonably homogeneous colors and albedos (see methods in 
Section 2). The families that passed the preliminary tests 
were scrutinized further to establish the best cutoff value 
dcut (the values of dcut were informed from past publica¬ 
tions), apply the V-shape criterion to identify large interlop¬ 
ers (section 4), and estimate tage from the envelope-based 
method (Section 5). 

We made sure that each identifled group with dcut was 
complete (i.e., not missing any extension that would make 
sense dynamically), and not part of a larger structure. In 
some cases, a cut in proper element space was applied to 
avoid the problem of chaining. In several cases, following 
Masiero et al. (2013), we used an albedo cutoff, py < 0.15, 
to identify dark families in the inner and central belts (the 
(84) Klio, (313) Chaldaea, (2262) Mitidika families). Ei- 
nally, we removed all identifled families and searched for 
signiflcant groups in the leftover population. This led to 
the identiflcation of several additional families (e.g., the 
(589) Croatia, (926) Imhilde, (1332) Marconia families in 
the outer belt). 

A difficulty we had to deal with is that different authors 
often listed the same family under different names. Eor ex¬ 
ample, Masiero et al. (2013) renamed the families after 
their largest fragments (e.g., the (158) Koronis family was 
referred to as the (208) Lacrimosa family). While it is very 
useful to identify the largest fragment, this naming con¬ 
vention creates chaos because some well-established family 
names disappear (Koronis, Maria, Elora, Geflon, Dora, etc.) 
and new names are used to refer to the same groups (Lac¬ 
rimosa, Roma, Augusta, Gudiachvili, Zhongolovich, etc.). 
It is also not clear in many cases what the true largest mem¬ 
ber is ((208) Lacrimosa is only slightly larger than (158) 
Koronis, some large objects have peripheral orbits and can 
be interlopers, etc.). 

A simple solution to this problem is to think of the family 
names as of labels, used from historical and other reasons, 
which do not necessarily reveal what the largest or lowest- 
numbered asteroid is. Thus, the Koronis family would re¬ 
main the Koronis family, regardless of whether (158) Ko¬ 
ronis, (208) Lacrimosa or some other large asteroid are the 
actual family members, largest or not. 

Anticipating that some authors will still prefer to rename 
some families, here we assign each family a unique Eamily 
Identifler Number (EIN). The EIN goes from 101 to 399 for 
families identifled from the analytic proper elements: IXX 
for the inner belt (2-2.5 AU), 2XX for the middle belt (2.5- 
2.82 AU), and 3XX for the outer belt (2.82-3.7 AU). The 
families identifled from the synthetic proper elements were 
assigned the EIN from 401 to 999: 4XX and 7XX for the 
low and high inclination families in the inner belt, 5XX and 
8 XX for the low and high inclination families in the mid¬ 
dle belt, and 6 XX and 9XX for the low and high inclination 
families in the outer belt. The division between low and 
high inclinations was taken at sin ip = 0.3 (ip 17.5°) 
which coincides with the highest inclinations for which the 
analytic proper elements are available. The EINs of the low- 
inclination families are ’aligned’ such that 101 and 401,102 
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Fig. 9.— The orbital location of notable families from Table 2. A triangle is placed at the orbit of an asteroid after which the family is 
named. The label near the triangle shows the designation number of that asteroid. 


and 402, 201 and 501, etc., are the same families identified 
from two different datasets. The families among Hungarias 
with a < 2 AU, Hildas with a 3.8 AU, and Jupiter Tro¬ 
jans with a b.2 AU are given the FIN between 001 and 
099. 

The main idea behind this notation is that even if a fam¬ 
ily is eventually renamed, its FIN will remain the same. 
Therefore, it should suffice to list the FIN of the renamed 
family in a new publication to make it clear how the new 
name relates to the designation(s) used in the past. For ex¬ 
ample, the Koronis family is given the FIN equal to 605 
(synthetic proper elements, outer belt, low i, 5th notable 
family), and will remain 605 even if the name Lacrimosa 
family will eventually be adopted. Since, however, it can 


be easily checked that the Koronis and Lacrimosa families 
have the same FIN (assuming that authors list the FIN in 
their publications), it will be crystal clear that we deal the 
same family. 

A potential difficulty with this scheme would arise if 
families that were previously assigned a FIN are often 
found to be nonexistent (say, because they are later found 
statistically insignificant or part of a large structure that was 
not evident in the original data). We suggest that such fam¬ 
ilies are deleted and the associated FINs are left void (i.e., 
not assigned to any family) such that there is no source for 
confusion. Finally, to avoid having many void FINs and not 
to run out of the FINs with the three digit scheme described 
above, the FIN should be given only to families that are 
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reasonably secure. Other, less secure families can be tenta¬ 
tively listed as family candidates and can wait for their FIN 
until the situation becomes clear with the new data. If the 
number of families in the inner, middle and/or outer belts 
exceeds 99, the 3-digit FIN scheme described here can be 
easily modified to a 4-digit scheme by adding 0 at the sec¬ 
ond decimal digit, such that, e.g., FIN=401 (the Vesta fam¬ 
ily) becomes FIN=4001. This change can be implemented 
if the need arises. 

7.2 FAMILY LISTS 

There are 122 notable asteroid families reported in Ta¬ 
ble 2 and 19 family candidates in Note 11. The vast ma¬ 
jority of families listed in Table 2 have been reported in 
many past publications and are clear beyond doubt (Vesta, 
Massalia, Eunomia, Gefion, Koronis, Eos, Themis, etc.). 
Many additional families, reported in only one or two pre¬ 
vious publications, are also clear beyond doubt with the 
new data. Eleven families are reported here for the first 
time ((329) Svea and (108138) 2001 GBll in the inner belt, 
(89) Julia, (369) Aeria, (1484) Postrema, (2782) Leonidas 
and (3152) Jones in the central belt, (589) Croatia, (926) 
Imhilde, (1332) Marconia and (106302) 2000 UJ87 in the 
outer belt). There are 20 notable families in the inner belt 
(the Nysa-Polana complex is counted as 3 families here), 47 
in the central belt, 46 in the outer belt, 2 among Hildas in 
the 3:2 resonance with Jupiter, 1 in Hungarias, and 6 among 
Jupiter Trojans. Of these, seventeen have high orbital incli¬ 
nation (ip > 17.5°). Figure 9 shows the orbital location of 
the notable families in the main belt. 

The lists of members of the notable families can be 
downloaded from the PDS node (Note 10). Each list con¬ 
tains the following information: (1) the asteroid number, 
(2) ap, (3) ep, (4) sin ip, and (5) the absolute magnitude H 
from the Minor Planet Center. Also, for families that have 
the V-shape envelope in (ap,i7), we fit Cq to this shape 
(Section 4, column 7 in Table 2) and report Cj / Cq for each 
family member in column 6 of the PDS files. The average 
albedo of each family was obtained from the WISE, and is 
reported in column 9 of Table 2 (py). The taxonomic type 
of families, reported in column 8 of Table 2, was taken from 
the previous taxonomic classification of families (Cellino et 
al. 2002) or was deduced from the SDSS colors. 

Columns 5 and 6 of Table 2 report the estimated diam¬ 
eter of the largest member, UpM, and diameter of a sphere 
with volume equivalent to that of all fragments, Dfrag. Dlm 
was obtained from AKARI, if available, or from WISE, if 
available, or was estimated from H and average py. The 
largest member and suspected interlopers with | Cj /Co | > 2 
were excluded in the estimate of Dfrag. The comparison 
of Dlm and Dfrag helps to establish whether a particular 
breakup event was catastrophic (Dfrag > Dlm) or crater¬ 
ing (Dfrag < Dlm), but note that this interpretation may 
depend on sometimes uncertain membership of the largest 
family objects. Also, the diameter of the parent body of 
each family can be estimated as Dpb = (^lm + 


but note that this estimate ignores the contribution of small 
(unobserved) fragments. 

7.3 COMPARISON WITH PREVIOUS DATASETS 

The family synthesis presented here is consistent with 
the results reported in Nesvorny (2012), Broz et al. (2013) 
and Carruba et al. (2013a). Eor example, all families re¬ 
ported in Nesvorny (2012) were found to be real here (ex¬ 
cept the (46) Hestia family that was moved to candidates; 
see Note 11). Nesvorny (2012), however, used very conser¬ 
vative criteria for the statistical significance of a family, and 
reported only 76 families (or 78 if the Nysa-Polana com¬ 
plex is counted, as it should, as three families). Using the 
2014 catalog of proper elements, albedo information from 
Masiero et al. (2013), and validating several new families 
from Milani et al. (2014), here we collect 44 families that 
did not appear in Nesvorny (2012). Almost all families re¬ 
ported in Broz et al. (2013) also appear here (a notable ex¬ 
ception is a large group surrounding (1044) Teutonia, which 
we believe is not a real family; see Note 11), but many new 
cases were added. The correspondence with Carruba et al. 
(2013a) is also good. 

Parker et al. (2008) used Nmin = 100 and therefore 
missed many small families which did not have more than 
100 members in the 2008 catalog. Also, given that they 
used a subset of asteroids with the SDSS colors, even rel¬ 
atively large families were unnoticed in this work (e.g. the 
(752) Sulamitis family in the inner belt now has 303 mem¬ 
bers). The high-i families were not reported in Parker et al. 
(2008), because they only used the analytic proper elements 
which are not available for the high-/ orbits. The strength 
of Parker et al.’s identification scheme was the reliability. 
Indeed, of all families reported in Parker et al. (2008) only 
the (1044) Teutonia, (1296) Andree and (2007) McCuskey 
(part of the Nysa-Polana complex) are not included among 
the notable families here (Parker et al.’s (110) Lydia family 
appears here as the (363) Padua family). 

Masiero et al. (2013; hereafter Ml3) reported 28 new 
cases and found that 24 old families were lost when com¬ 
pared to the family lists in Nesvorny (2012). Most families 
not listed in Ml3 are well-defined families such as Karin, 
Beagle, Datura, Emilkowalski, Lucascavin, etc. These fam¬ 
ilies were not listed in Ml3, because they overlap with 
larger families (and were included in their membership lists 
in Ml3) or because they only have a few known members 
(i.e., fall below A/min used in Ml3). On the other hand, 
we verified that many new cases reported in Masiero et al. 
(2013) are genuine new families that can be conveniently 
found with the albedo cutoff (e.g., the (84) Klio, (144) Vi- 
bilia, (313) Chaldaea (= (1715) Salli in M13), (322) Phaeo, 
(623) Chimarea, (816) Juliana, and (1668) Hanna families). 
These families were included here. In some cases, we found 
that M13’s new family barely stands out from the back¬ 
ground and it thus seems uncertain. To stay on the safe side, 
we therefore report these cases as the candidate families in 
Note 11. These families may be real, but their statistical 
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significance needs to be carefully tested with the new data. 

Finally, we compare the family synthesis with Milani et 
al. (2014; hereafter Ml4), who used the newest catalogs 
from all previous works discussed here. They identified 
many new families which are certain beyond doubt (e.g., 
the (96) Aegle, (618) Elfriede, (2344) Xizang, (3438) Inar- 
radas, (3811) Karma, (7468) Anfimov, (53546) 2000 BY6; 
the (96) Aegle and (618) Elfriede families were also re¬ 
ported in Ml3). These cases highlight the strength of the 
Ml4 identification scheme and are included in the family 
synthesis in Table 2. Some smaller families located inside 
bigger families (e.g., the (832) Karin family families in the 
(158) Koronis family, (656) Beagle in (24) Themis) were 
not reported in Ml4. In addition, several families were 
not reported, presumably because the QRL was set too low 
to detect them. This happens, most notably, in the 2.82- 
2.96 AU region (i.e., between the 5:2 and 7:3 resonances), 
where the number density of asteroids is relatively low (see 
Figures 1 and 9). The notable cases in this region in¬ 
clude the (81) Terpsichore, (709) Fringilia, (5567) Durisen, 
(5614) Yakovlev, (7481) San Marcello and (10811) Lau 
families. A possible solution to this issue would be, in terms 
of the M14 identification scheme, to use a separate QRL 
level for the 2.82-2.96 AU region. Also, several families 
were split in M14 into several parts. This affects the follow¬ 
ing families: (8) Flora (split in 4 parts), (31) Euphrosyne (3 
parts), (221) Eos (5 parts), (702) Alauda (4 parts), (1400) 
Tirela (8 parts), (2085) Henan (4 parts), (4203) Brucato (4 
parts) (here we only list families that were split to three or 
more parts in Ml4). 

8. CONCLUSIONS 

It is clear from several different arguments that the list 
of known families must be largely incomplete. For exam¬ 
ple, most families with an estimated parent body size below 
100-200 km are found to have ages tage ^ 1 g-y* (^-g-^ 
Broz et al., 2013). In contrast, the rate of impacts in the 
main belt, and therefore the rate of family-forming events, 
should have been roughly unchanging with time over the 
past :^3.5 g.y. (and probably raised quite a bit toward the 
earliest epochs). So, there must be many missing families 
with the formation ages tage > 1 g-Y- (Note 12). These 
families are difficult to spot today, probably because they 
have been dispersed by dynamical processes, lost members 
by collisional grinding, and therefore do not stand out suf¬ 
ficiently well above the dense main-belt background (they 
are now part of the background). 

The significant incompleteness of known families is also 
indicated by the extrapolation of the number of families de¬ 
tected in the 2.82-2.96 AU zone to the whole main belt. 
As we hinted on at the end of the last section, the 2.82- 
2.96 AU zone is sparsely populated such that asteroid fam¬ 
ilies can be more easily detected there. Nearly 20 families 
with ip < 17.5° were identified in this region. In com¬ 
parison, the 2.96-3.3 AU zone, where also :^20 families 
with ip < 17.5° can be found, is about twice as wide as 


the 2.82-2.96 AU zone and contains about twice as many 
large asteroids. A straightforward conclusion that can be 
inferred from this comparison, assuming everything else be¬ 
ing equal, is that the families in the 2.96-3.3 AU zone are 
(at least) a factor of ^2 incomplete. A similar argument 
applies to the inner and central belts. 

This is actually good news for future generations of plan¬ 
etary scientists, because this field is open for new discover¬ 
ies. Figuring out how to find the missing asteroid families 
with fage > 1 g-Y- will not be easy. One way forward 
would be to improve our capability to model the dynam¬ 
ical evolution of main belt asteroids over g.y. timescales, 
such that we can rewind the clock and track fragments back 
to their original orbits. The modeling of the Yarkovsky ef¬ 
fect could be improved, for example, if we knew the spin 
states, densities, conductivity, etc., of small main-belt aster¬ 
oids on individual basis. Another approach would consist 
in identifying families based on the physical properties of 
their members. While this method is already in use, mainly 
thanks to data from the SDSS and WISE, we anticipate that 
it can be pushed much further, say, with automated spec¬ 
troscopic surveys, or, in more distant future, with routine 
sampling missions. 

9. NOTES 

Note 1 - Alternatively, one can use the frequencies n, g and s (Carruba 
and Michtchenko, 2007, 2009), where n is the mean orbital frequency and 
g and s are the (constant) precession frequencies of the proper perihelion 
longitude tup and the proper nodal longitude tip, respectively. The use 
of frequencies, instead of the proper elements, can be helpful for asteroid 
families near or inside the secular orbital resonances (e.g. the Tina family 
in the uq resonance; Carruba and Morbidelli, 2011). 

Note 2 - Additional methods were developed and/or adapted for spe¬ 
cific populations of asteroids such as the ones on the high-inclination and 
high-eccentricity orbits (Lemaitre and Morbidelli, 1994) or in orbital res¬ 
onances (Morbidelli, 1993; Milani, 1993; Beauge and Roig, 2001; Broz 
and Vokrouhlicky, 2008; Broz and Rozehnal, 2011). 

Note 3 - http : //hamilton . dm. unipi . it/astdys/ 

Note 4 - See chapter Bendjoya and Zappala (2002) in the Asteroids III 
book for a discussion of other clustering algorithms such as the Wavelet 
Analysis Method (WAM) and D-criterion. The WAM was shown to pro¬ 
duce results that are in a good agreement with those obtained from the 
HCM (Zappala et al., 1994). The D-criterion was originally developed to 
identify meteorite streams (Southwork and Hawkins, 1963). These meth¬ 
ods have not been used for the classification of asteroid families in the past 
decade, and we therefore do not discuss them here. 

Note 5 - The SDSS measured flux densities in five bands with effective 
wavelengths 354, 476, 628, 769, and 925 nm. The WISE mission mea¬ 
sured fluxes in four wavelengths (3.4, 4.6, 12 and 22 gm), and combined 
the measurements with a thermal model to calculate albedos (pv) and 
diameters (D). The latest public releases of these catalogs include color 
or albedo data for over 100,000 main belt asteroids with known orbits, of 
which about 25,000 have both the color and albedo measurements. The 
catalogs are available at 

WWW. sdss . org/dr 6/product s/value-added/index . html 
and irsa.ipac.caltech.edu/Missions/wise.html. 

Note 6 - Most but not all asteroid families are physically homogeneous. 
The Eos family has the highest diversity of taxonomic classes of any 
known family (e.g., Mothe-Diniz et al., 2008). This diversity has led to 
the suggestion that the Eos parent body was partially differentiated. It can 
also be the source of carbonaceous chondrites (Clark et al., 2009). The 
Eunomia family may be another case of a relatively heterogeneous fam¬ 
ily (e.g., Nathues et al., 2005). See Weiss and Elkins-Tanton (2013) for a 
review. 
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Note 7 - DeirOro et al. (2004) attempted to model observed family shapes 
by Gaussian ellipsoids. The distribution of |/| obtained in this work was 
strongly peaked near 7r/2, while a more uniform distribution between 0 
and TT would be expected if different breakups occurred at random orbital 
phases. This result was attributed to the Yarkovsky effect. 

Note 8 - Nesvorny et al. (2008a) found evidence for a large population 
of V-type asteroids with slightly lower orbital inclinations (ip = 3-4°) 
than the Vesta family (zp ~ 5°). Because these asteroids could not have 
dynamically evolved from the Vesta family region to their present orbits in 
~1 g.y., they are presumably fragments excavated from (4) Vesta’s basaltic 
crust by an earlier impact. 

Note 9 - Other asteroid fa mili es whose long-term dynamics has been stud¬ 
ied in detail, listed here in the alphabetical order, are the Adeona family 
(affected in ep and ip by the 8:3 resonance at 2.705 AU, Carruba et al., 
2003), Agnia family (inside the zi resonance; Vokrouhlicky et al., 2006c), 
Astrid family (near the border of the 5:2 resonance; Vokrouhlicky et al., 
2006b), Eunomia family (Carruba et al., 2007b), Euphrosyne family (lo¬ 
cated in a region with many resonances, including g — ge = 0, near the 
inner border of the 2:1 resonance; Carruba et al., 2014), Erigone family 
(cut in the middle by the Z 2 resonance; Vokrouhlicky et al., 2006b), Gefion 
family (affected by various resonances near 2.75 AU, Carruba et al., 2003, 
Nesvorny et al., 2009), Hilda and Schubart families in the 3:2 resonance 
with Jupiter (Broz and Vokrouhlicky, 2008), Hungaria family (perturbed 
by 2g — g^— ge = 0 and other secular resonances below 1.93 AU; Warner 
et al. (2009), Milan! et al. (2010), also see Galiazzo et al. (2013, 2014) for 
the contribution of Hungarias to the E-type NEAs and Cuk et al. (2014) for 
their suggested relation to the aubrite meteorites), Hygiea family (Carruba, 
2013; Carruba et al., 2014), Massalia family (the part with ap > 2.42 AU 
disturbed by the 1:2 resonance with Mars, Vokrouhlicky et al., 2006b), 
Merxia family (spread by the 3J-1S-1 three-body resonance at up = 2.75 
AU; Vokrouhlicky et al., 2006b), Padua family (Carruba, 2009a), Pallas 
family (Carruba et al., 2011), Phocaea family (Carruba, 2009b), Sylvia 
family ((87) Sylvia has two satellites, possibly related to the impact that 
produced the Sylvia family, Vokrouhlicky et al., 2010), and Tina family 
(Carruba and Morbidelli, 2011). 

Note 10-http ://sbn.psi.edu/pds/resource/nesvornyf am 

Note 11 - The candidate families are: (929) Algunde, (1296) Andree, 
(1646) Rosseland, (1942) Jablunka, (2007) McCuskey, (2409) Chapman, 
(4689) Donn, (6246) Komurotoru and (13698) 1998 KF35 in the inner belt, 
(46) Hestia, (539) Palmina, (300163) P/2006 VW 139, (3567) Alvema, 
(7744) 1986 QAl in the central belt, and (260) Huberta, (928) Hilrun, 
(2621) Goto, (1113) Katja, (8737) Takehiro in the outer belt. We ten¬ 
tatively moved the (46) Hestia family, previously known as FIN 503, to 
the family candidate status, because this group is not convincing with the 
present data. The previously reported groups around (5) Astraea, (1044) 
Teutonia, (3110) Wagman, (4945) Ikenozenni, (7744) 1986 QAl, (8905) 
Bankakuko, (25315) 1999 AZ8, (28804) 2000 HC81 seem to align with 
the 2 ;i = g s — ge — sq = 0 resonance and are probably an artifact of 
the HCM chaining. 

Note 12 - The list of known families corresponding to parent bodies with 
D > 200 km, on the other hand, is probably reasonably complete, because 
the estimated ages of these families appear to be randomly distributed over 
4 g.y. (Broz et al., 2013). These largest families can therefore be used to 
constrain the collisional history of the asteroid belt (see chapter by Bot- 
tke et al. in this volume). 
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Table 2: Notable asteroid families. The columns are the: (1) Family Identification Number (FIN), (2) family name, (3) cutoff distance (dcuu asterisk denotes dcut 
used on a subset of asteroids with py < 0.15), (4) number of family members identified with dcut. (5) largest member(s) in the family (either the number designation 
of the largest member(s), in parenthesis, if different from the asteroid after which the family is named, or the estimated diameter of the largest member, T^lm), ( 6 ) 
diameter of a sphere with volume equivalent to that of all fragments (I^frag, (7) Co parameter defined in Section 4, (8) taxonomic type, (9) mean geometric albedo 
from WISE (py), and (10) various references and notes. We do not report t^ge here but note that t^ge can be estimated from Co given in column 7 and Eq. (2). See 
http : // sirrah . tro ja . mf f . cuni . cz/ ~mira/mp/f ams/ for additional information on asteroid families (Broz et al., 2013). 


FIN 

Family 

Name 

dcut 

(ms-i) 

#of 

mem. 

-Dlm 

(km) 

-t^frag 

(km) 

1 

o 

T-1 

Tax. 

Type 

py 

References and Notes 

001 

153 Hilda 

130 

409 

Hungarias, Hildas and Jupiter Trojans 

164 - - C 0.04 

Broz et al. (2011) 

002 

1911 Schubart 

60 

352 

80 

91 

- 

c 

0.03 

Broz and Vokrouhlicky (2008) 

003 

434 Hungaria 

100 

2965 

10 

24 

0.3 zb 0.1 

E 

0.35 

Warner et al. (2010), Milan! et al. (2011) 

004 

624 Hector 

50 

12 

231 

- 

- 

- 

- 

satellite, Marchis et al. (2014), Rozehnal and Broz (2014) 

005 

3548 Eurybates 

50 

218 

68 

87 

- 

CP 

0.06 

Roig et al. (2008), Broz and Rozehnal (2011) 

006 

9799 1996 RJ 

60 

7 

72 

26 

- 

- 

0.06 

Rozehnal and Broz (2014) 

007 

James Bond 

oo 

1 

(himself) 

- 

- 

ASP 

variable 

Campbell et al. (1995) 

008 

20961 Arkesilaos 

50 

37 

- 

- 

- 

- 

- 

Rozehnal and Broz (2014) 

009 

4709 Ennomos 

100 

30 

(1867,4709) 

- 

- 

- 

0.06 

Rozehnal and Broz (2014) 

010 

247341 2001 UV209 

100 

13 

- 

- 

- 

- 

0.09 

Rozehnal and Broz (2014) 

401 

4 Vesta 

50 

15252 

Inner Main Belt, 2.0 < a < 2.5 AU, i 

525 50 1.5 zb 0.5 V 

< 17.5° 

0.35 

source of HEDs, two overlaping families? 

402 

8 Elora 

60 

13786 

(8,254) 

- 

2.5 zb 0.5 

S 

0.30 

dispersed, source of LL NEAs, Dykhuis et al. (2014) 

403 

298 Baptistina 

45 

2500 

21 

- 

0.25 zb 0.05 

X 

0.16 

related to K/T impact? Bottke et al. (2007) 

404 

20 Massalia 

55 

6424 

132 

27 

0.25 zb 0.05 

S 

0.22 

Vokrouhlicky et al. (2006b) 

405 

44 Nysa-Polana 

50 

19073 

(135,142,495) 

- 

1.0 zb 0.5 

SEC 

0.28/0.06 

Walsh et al. (2013), Dykhuis and Greenberg (2015) 

406 

163 Erigone 

50 

1776 

72 

46 

0.2 zb 0.05 

CX 

0.06 

Vokrouhlicky et al. (2006b) 

407 

302 Clarissa 

55 

179 

34 

15 

0.05 zb 0.01 

X 

0.05 

compact with ears, cratering 

408 

752 Sulamitis 

55 

303 

61 

35 

0.3 zb 0.1 

C 

0.04 


409 

1892 Lucienne 

100 

142 

11 

11 

0.15 zb 0.05 

S 

0.22 


410 

27 Euterpe 

65 

474 

no 

16 

0.50 zb 0.25 

S 

0.26 


411 

1270 Datura 

10 

6 

8 

3 

- 

S 

0.21 

Nesvorny et al. (2006c) 

412 

21509 Lucascavin 

10 

3 

- 

- 

- 

S 

- 

Nesvorny and Vokrouhlicky (2006) 

413 

84 Klio 

130* 

330 

78 

33 

0.75 zb 0.25 

c 

0.07 

interloper 12?, Masiero et al. (2013) 

414 

623 Chimaera 

120 

108 

43 

21 

0.3 zb 0.1 

CX 

0.06 

Masiero et al. (2013) 

415 

313 Chaldaea 

130* 

132 

(313,1715) 

- 

1.0 zb 0.5 

c 

0.07 

1715 in Masiero et al. (2013) 

416 

329 Svea 

150 

48 

70 

21 

0.3 zb 0.1 

CX 

0.06 

new, near 3:1 

417 

108138 2001 GBll 

20 

9 

- 

- 

- 

- 

- 

new, compact 

701 

25 Phocaea 

150 

1989 

Inner Main Belt, 2.0 < a < 2.5 AU, 
(25,587) - 2.0 zb 1.0 S 

i > 17° 

0.22 

Carruba (2009b), Carruba et al. (2010) 





Table 2: Continued. 


FIN 

Family 

Name 

dcut 

(ms-i) 

# of 

mem. 

D-lm 

(km) 

-(^frag 

(km) 

1 

o 

T- 1 

Tax. 

Type 

pv 

References and Notes 

501 

3 Juno 

55 

1684 

231 

Central Main Belt, 2.5 < a < 2.S2AU, 

25 0.5 zb 0.2 S 0.25 

i < 17.5° 

cratering, relation to H chondrites? 

502 

15 Eunomia 

50 

5670 

256 

100 

2.0 zb 0.7 

s 

0.19 

continues beyond 5:2? 

503 

- 

- 

- 

- 

- 

- 

- 

- 

46 Hestia moved to candidates 

504 

128 Nemesis 

50 

1302 

178 

50 

0.25 zb 0.05 

c 

0.05 

3827 in Milani et al. (2014), 125 in Cellino et al. (2002) 

505 

145 Adeona 

50 

2236 

141 

78 

0.7 zb 0.3 

c 

0.07 


506 

170 Maria 

60 

2940 

(472,170) 

- 

2.0 zb 1.0 

s 

0.25 

(472) Roma in Masiero et al. (2013) 

507 

363 Padua 

45 

1087 

91 

48 

0.5 zb 0.2 

X 

0.10 

Carruba (2009a), also known as the (110) Lydia family 

508 

396 Aeolia 

20 

296 

46 

13 

0.075 zb 0.025 

X 

0.17 

compact, young? 

509 

410 Chloris 

80 

424 

107 

56 

0.75 zb 0.25 

c 

0.06 

eroded 

510 

569 Misa 

50 

702 

65 

57 

0.5 zb 0.2 

c 

0.03 

V-shaped subfamily inside 

511 

606 Brangane 

55 

195 

36 

18 

0.04 zb 0.01 

s 

0.10 

compact, 606 offset, interloper? 

512 

668 Dora 

45 

1259 

(1734,668) 

- 

- 

c 

0.05 

668 offset, 1734 in Masiero et al. (2013), V-shaped subfamily 

513 

808 Merxia 

55 

1215 

34 

28 

0.3 zb 0.1 

s 

0.23 

Vokrouhlicky et al. (2006b) 

514 

847 Agnia 

30 

2125 

(847,3395) 

- 

0.15 zb 0.05 

s 

0.18 

resonance, Vokrouhlicky et al. (2006c) 

515 

1128 Astrid 

60 

489 

42 

29 

0.12 zb 0.02 

c 

0.08 

Vokrouhlicky et al. (2006b) 

516 

1272 Gefion 

50 

2547 

(2595,1272) 

- 

0.8 zb 0.3 

s 

0.20 

source of L chondrites? Nesvorny et al. (2009), also known as 93 and 2595 

517 

3815 Konig 

55 

354 

22 

34 

0.06 zb 0.03 

cx 

0.04 

compact, young? Nesvorny et al. (2003), 342 and 1639 offset 

518 

1644 Rafita 

70 

1295 

(1658,1587) 

- 

0.5 zb 0.2 

s 

0.25 

1644 probably interloper 

519 

1726 Hoffmeister 

45 

1819 

(272,1726) 

- 

0.20 zb 0.05 

CP 

0.04 

(272) Antonia in Masiero et al. (2013), but 272 offset 

520 

4652 lannini 

25 

150 

5 

10 

- 

s 

0.32 

1547 offset, compact, Nesvorny et al. (2003) 

521 

7353 Kazuya 

50 

44 

11 

10 

- 

s 

0.21 

small clump 

522 

173 Ino 

50 

463 

161 

21 

0.5 zb 0.2 

s 

0.24 

also known as 18466, large and dark 173 is probably interloper, ears? 

523 

14627 Emilkowalski 

10 

4 

7 

3 

- 

s 

0.20 

Nesvorny and Vokrouhlicky (2006) 

524 

16598 1992 YC2 

10 

3 

- 

- 

- 

s 

- 

Nesvorny and Vokrouhlicky (2006) 

525 

2384 Schulhof 

10 

6 

12 

4 

- 

s 

0.27 

Vokrouhlicky and Nesvorny (2011) 

526 

53546 2000 BY 6 

40 

58 

8 

18 

- 

c 

0.06 

Milani et al. (2014) 

527 

5438 Lorre 

10 

2 

30 

- 

- 

c 

0.05 

Novakovic et al. (2012) 

528 

2782 Leonidas 

50 

135 

(4793,2782) 

- 

- 

cx 

0.07 

new, related to 144? 

529 

144 Vibilia 

100* 

180 

142 

- 

- 

c 

0.06 

Masiero et al. (2013), PDS list identical to 2782 

530 

322 Phaeo 

100* 

146 

72 

31 

0.3 zb 0.1 

X 

0.06 

Cellino et al. (2002), joins (2669) Shostakovich 

531 

2262 Mitidika 

100* 

653 

(404,5079) 

- 

- 

c 

0.06 

dispersed, 404 offset, 2262 has pv = 0.21 

532 

2085 Henan 

50 

1872 

18 

32 

0.75 zb 0.25 

L 

0.20 

2085 offset in zp, 4 families in Milani et al. (2014) 

533 

1668 Hanna 

60 

280 

22 

32 

0.2 zb 0.1 

cx 

0.05 

Masiero et al. (2013) 

534 

3811 Karma 

60 

124 

26 

24 

0.25 zb 0.05 

cx 

0.05 

Milani et al. (2014) 
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FIN 

Family 

Name 

dcut 

(ms-i) 

#of 

mem. 

-Dlm 

(km) 

■^frag 

(km) 

o 

1 

Tax. 

Type 

pv 

References and Notes 

535 

2732 Witt 

45 

1816 

11 

Central Main Belt, 2.5 

25 0.75 zb 0.25 

< a < 

S 

2.32 AU, i < 17.5° 

0.26 relation to the Charts family beyond 5:2? 10955 and 19466 in Milani et al. 

536 

2344 Xizang 

65 

275 

17 

20 

0.3 zb 0.1 

- 

0.12 

Milani et al. (2014), includes 396 

537 

729 Watsonia 

130 

99 

52 

38 

- 

L 

0.13 

Cellino et al. (2014) 

538 

3152 Jones 

40 

22 

37 

11 

- 

T 

0.05 

new, compact, diagonal in (ap, ep) 

539 

369 Aeria 

90 

272 

75 

17 

0.3 zb 0.1 

X 

0.17 

new, part above 2.6778 AU down in zp 

540 

89 Julia 

70 

33 

147 

6 

- 

S 

0.19 

new, compact 

541 

1484 Postrema 

100 

108 

47 

- 

- 

cx 

0.05 

new, large 387,547,599? 

801 

2 Pallas 

350 

128 

513 

Central Main Belt, 2.5 

40 

< a < 

B 

2.32 AU, i > 17.5° 

0.16 Carruba et al. (2010, 2012), part beyond 5:2 

802 

148 Gallia 

200 

182 

81 

19 

0.5 zb 0.1 

S 

0.17 

large interlopers 

803 

480 Hansa 

200 

1094 

56 

62 

- 

S 

0.26 

2 families in Milani et al. (2014) 

804 

686 Gersuind 

120 

415 

49 

36 

- 

S 

0.15 

2 families in Milani et al. (2014) 

805 

945 Barcelona 

150 

306 

27 

19 

0.25 zb 0.05 

S 

0.25 

2 families in Milani et al. (2014) 

806 

1222 Tina 

200 

96 

26 

10 

0.10 zb 0.05 

X 

0.34 

in the ^ — ^6 = 0 resonance, Carruba and Morbidelli (2011) 

807 

4203 Brucato 

200 

342 

18 

44 

0.5 zb 0.2 

cx 

0.06 

1263 interloper? Carruba (2010), 4 families in Milani et al. (2014) 

601 

10 Hygiea 

60 

4854 

428 

Outer Main Belt, 2.82 < a < 3.7 AU, i < 17° 

- CB 0.06 Carruba et al. (2014) 

602 

24 Themis 

60 

4782 

177 

230 

2.5 zb 1.0 

c 

0.07 

includes 656 Beagle, Nesvorny et al. (2008b) 

603 

87 Sylvia 

130 

255 

263 

- 

- 

X 

0.05 

Vokrouhlicky et al. (2010) 

604 

137 Meliboea 

85 

444 

(511,137) 

- 

- 

c 

0.05 

(511) Davida in Masiero et al. (2013) 

605 

158 Koronis 

45 

5949 

(208,158,462) 

- 

2.0 zb 1.0 

s 

0.15 

(208) Lacrimosa in Masiero et al. (2013) 

606 

221 Eos 

45 

9789 

(221,579,639) 

- 

1.5 zb 0.5 

K 

0.13 

Vokrouhlicky et al. (2006a), Broz and Morbidelli (2013) 

607 

283 Emma 

40 

76 

122 

56 

0.3 zb 0.1 

C 

0.05 

affected by the zi resonance? 

608 

293 Brasilia 

50 

579 

(3985) 

- 

0.20 zb 0.05 

X 

0.18 

293 interloper?, also known as (1521) Sejnajoki, Nesvorny et al. (2003) 

609 

490 Veritas 

30 

1294 

113 

78 

0.2 zb 0.1 

CPD 

0.07 

see Section 3 

610 

832 Karin 

10 

541 

17 

16 

0.03 zb 0.01 

S 

0.21 

see Section 3, Harris et al. (2009) 

611 

845 Naema 

35 

301 

61 

37 

0.20 zb 0.05 

c 

0.08 


612 

1400 Tirela 

50 

1395 

(1040,1400) 

- 

0.75 zb 0.25 

s 

0.07 

8 families in Milani et al. (2014) 

613 

3556 Lixiaohua 

45 

756 

(3330,3556) 

- 

0.25 zb 0.05 

cx 

0.04 

3330 offset, Novakovic et al. (2010) 

614 

9506 Telramund 

45 

468 

(179,9506) 

- 

- 

s 

0.22 

(179) Klytaemnestra in Masiero et al. and Milani et al. 

615 

18405 EY12 

50 

104 

9 

14 

0.08 zb 0.03 

cx 

0.17 
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FIN 

Family 

Name 

dcut 

(ms-i) 

#of 

mem. 

-Dlm 

(km) 

■^frag 

(km) 

Co 

(10-^ AU) 

Tax. 

Type 

pv 

References and Notes 

616 

627 Charis 

80 

808 

Outer Main Belt, 2.82 < a < 3.7 AU, i < 17.5° 

50 45 - C 0.08 

16286 in Milani et al. (2014), related to Witt family? 

617 

778 Theobalda 

60 

376 

66 

50 

- 

cx 

0.06 

778 offset, Novakovic et al. (2010) 

618 

1189 Terentia 

60 

79 

63 

18 

0.13 zb 0.03 

c 

0.07 


619 

10811 Lau 

120 

56 

8 

6 

0.075 zb 0.025 

s 

0.27 

36824 interloper? 

620 

656 Beagle 

25 

148 

54 

28 

0.07 zb 0.03 

c 

0.09 

656 and 90 offset, Nesvorny et al. (2008b) 

621 

158 Koronis(2) 

8 

246 

34 

13 

0.010 zb 0.005 

s 

0.14 

young, Molnar and Haegert (2009) 

622 

81 Terpsichore 

120 

138 

123 

27 

0.50 zb 0.25 

c 

0.05 


623 

709 Fringilla 

150 

134 

96 

55 

- 

X 

0.05 

large 1191 

624 

5567 Durisen 

100 

27 

34 

23 

- 

X 

0.04 


625 

5614 Yakovlev 

120 

67 

13 

23 

0.15 zb 0.05 

c 

0.05 


626 

7481 San Marcello 

100 

144 

(3978,7489) 

- 

- 

X 

0.19 

also known as 12573 

627 

15454 1998 YB3 

80 

38 

(3156,15454) 

- 

0.10 zb 0.05 

cx 

0.05 


628 

15477 1999 CGI 

95 

248 

- 

- 

- 

s 

0.10 


629 

36256 1999 XT17 

70 

58 

(15610,36256) 

- 

0.25 zb 0.05 

s 

0.21 

several large bodies 

630 

96 Aegle 

50 

99 

165 

38 

- 

cx 

0.07 

Masiero et al. and Milani et al. 

631 

375 Ursula 

70 

1466 

(1306,375) 

- 

- 

cx 

0.06 

Masiero et al. and Milani et al. 

632 

618 Elfriede 

40 

63 

122 

26 

- 

c 

0.05 

compact, recent? 

633 

918 Itha 

100 

54 

21 

35 

- 

s 

0.23 

dispersed, many sizable members 

634 

3438 Inarradas 

80 

38 

25 

33 

- 

cx 

0.07 

Milani et al. (2014) 

635 

7468 Anfimov 

60 

58 

10 

14 

- 

s 

0.16 

Milani et al. (2014) 

636 

1332 Marconia 

30 

34 

50 

16 

- 

cx 

0.05 

new 

637 

106302 2000 UJ87 

60 

64 

7 

15 

- 

cx 

0.05 

new, large 132999 

638 

589 Croatia 

45 

93 

92 

31 

0.5 zb 0.2 

X 

0.07 

new, 21885 in Milani et al. (2014) 

639 

926 Imhilde 

70 

43 

50 

18 

0.2 zb 0.1 

cx 

0.05 

new 

640 

P/2012 F5 (Gibbs) 

10 

8 

- 

- 

- 

- 

- 

Novakovic et al. (2014) 

641 

816 Juliana 

80 

76 

68 

39 

- 

cx 

0.05 

Masiero et al. (2013) 

901 

31 Euphrosyne 

120 

2035 

Outer Main Belt, 2.82 < a < 3.5 AU, i > 

276 130 - C 

17.5° 

0.06 

Carruba et al. (2014), 3 families in Milani et al. (2014) 

902 

702 Alauda 

120 

1294 

191 

- 

2.5 zb 1.0 

B 

0.07 

276 and 1901 offset, 4 families in Milani et al. (2014) 

903 

909 Ulla 

120 

26 

113 

28 

- 

X 

0.05 


904 

1303 Luther a 

50 

163 

87 

56 

- 

X 

0.04 

also known as (781) Kartivelia 

905 

780 Armenia 

50 

40 

98 

22 

- 

c 

0.05 

compact 



